DNA-TR-90-78 


Impact  of  Ionospheric  Scintillations 
on  Space-Based  Surveillance  Radars 


Lester  L.  DeRaad,  Jr. 
Brian  Lamb 
R  &  D  Associates 
P.O.  Box  92500 
Los  Angeies,  CA  90009 


June  1991 


Technical  Report 


DTIC 


ELECTE 

JULOSm 


B 


CONTRACT  No,  DNA  001-88-C-0046 


Approved  for  public  release: 
distribution  is  unlimited. 


91-03971 


Destroy  this  report  when  it  is  no  longer  needed.  Do  not 
return  to  sender. 


PLEASE  NOTIFY  THE  DEFENSE  NUCLEAR  AGENCY, 
ATTN:  CSTI,  6801  TELEGRAPH  ROAD,  ALEXANDRIA,  VA 
22310-3398,  IF  YOUR  ADDRESS  IS  INCORRECT,  IF  YOU 
WISH  IT  DELETED  FROM  THE  DISTRIBUTION  LIST,  OR 
IFTHE  ADDRESSEE  IS  NO  LONGER  EMPLOYED  BY  YOUR 
ORGANIZATION. 


CUT  HERE  AND  RETURN 


DISTRIBUTION  LIST  UPDATE 


This  mailer  is  provid  i  to  enable  DNA  to  maintain  current  distribution  lists  for  reports.  We  would 
appreciate  your  providing  the  requested  information 

□  Add  the  individual  listed  to  your  distribution  list. 

□  Delete  the  cited  organization/individual. 

□  Change  of  address. 

NAME: _ 

ORGANIZATION: _ 


OLD  ADDRESS  CURRENT  ADDRESS 


NOTE: 

Please  return  the  mailing  label  from 
the  document  so  that  any  additions, 
changes,  corrections  or  deletions 
can  be  made  more  easily. 


DNA  OR  OTHER  GOVERNMENT  CONTRACT  NUMBER- _ 

CERTIFICATION  OF  NEED-TO-KNOW  BY  GOVERNMENT  SPONSOR  (if  other  than  DNA): 

SPONSORING  ORGANIZATION: _ 

CONTRACTING  OFFICER  OR  REPRESENTATIVE: _ 


SIGNATURE 


Director 

Defense  Nuclear  Agency 
ATTN:  TITL 

Washington,  DC  20305-?. 


Director 

Defense  Nuclear  Agency 
ATTN:  TIT). 

Washington.  DC  20305-1000 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0168 


Public  reporting  burden  (or  this  collection  o*  information  ts  estimated  to  average  1  hour  per  response  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed  ard  completing  and  reviewing  the  collect'On  o*  mfon"  ioon  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  (his 
collection  of  information  including  suggest  jns  for  reducing  this  burden  to  Washington  Headqur^ile  s  Serviciis  Directorate  for  information  Operations  and  Reports.  1215  Jefferson 
Davis  H.ghway.  Suite  1204.  Arlington.  VA  22202-4302,  and  to  the  OlfiCv-  of  Management  3'id  ''udget.  Pape^ork  Reduction  Project  (0704-0138).  Washington.  DC  20503 


1  AGENCY  USE  ONLY  (Leave  blank) 


4.  riTLS  AND  SUBTITLE 


2  REPORT  DAT  c 

910601 


3.  REPORT  TYPE  AND  DATES  COVERED 

Technical  890 103-891231 


5  FUNDING  NUMBERS 


Impact  of  ionospheric  Sclntiliatlons  on  Space-Based 
Surveillance  Radars 


6.  AUTHOR'S) 

Les':er  L.  DeRaad.  Jr.,  and  Brian  Lamb 


C  -  DNA  001-88-C-0046 

PE  -  627 15H 

PR  -  RD 

TA  -  RC 

WU  -  DH041780 


7  PERFORMING  ORGANIZATION  NAME(S)  AND  ADORESS(ES) 

R  &  D  Associates 

P.O.  Box  92500 

Los  Angeles,  CA  90009 


8  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

RDA-TR- 161603-001 


10  SPONSORING/MONITOBING 
AGENCY  REPORT  NUMBER 


DNA-TR-90-78 


9  SPONSORING/MONi  'ORING  AGENCY  NAME(S)  AND  ADDneSS(ES) 

Defense  Nuclear  Agency 
6801  Telegraph  Road 
Alexandria,  VA  22310-3398 
RAAE/Ullrich 


11.  SUPPLEMENTARY  NOTES 

This  work  was  sponsored  by  the  Defense  Nuclear  Agency  under  RDT&E  RMC  Code  B4662D 
RD  RC  00021  DFPR  1910A  25904D. 


12a.  OISTRIBUTION/AVAILABILITY  STATEMENT 


12b  DISTRIBUTION  CODE 


Approved  for  public  release:  distribution  is  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

A  model  is  developed  to  analyze  propagation  effects  on  space-based  radars  (SBRs)  with 
pulse-Doppler  processing  and  with  a  displaced  phase  center  antenna  (DPCA)  for  rejection 
of  surface  clutter.  It  is  found  that  nuclear-disturbed  ionospheric  effects  can  significantly 
affect  radar  performance.  Depending  on  the  radar  design  parameters,  simlla  .  but  less  se¬ 
vere,  problems  may  cdso  occur  in  ambient  auroral  and  equatorial  ionospheres.  The  iono¬ 
spheric  effects  degrade  the  DPCA  clutter  rejection  and  target  Doppler  processing  funciions 
of  the  radar. 


Cl  in  ICr'T  TCDAXC 


Space-Based  Radar 
Propagation  Effects 


Rayleigh  Fading 
Scintillation  Effects 


15.  NUMBER  OF  PAGES 

70 


16  PRICE  CODE 


17  SECURITY  CLASSIFICAriON  18  SECURITY  CLASSIFICATION  19.  SECURI  Pv  CLASSIFICATION  20  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 


UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 


NSN  7540-280-5500 


1 


Standard  Form  298  (Rev  2-89) 
Prescribed  by  ANSI  Sla  239  1 8 
298  102 


secuqiTY  classification  of  this  page 

UNCLASSIFIED 


SUMMARY 
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It  is  found  that  ionospheric  propagation  effects  following  a  high-altitude  nuclear  explosion 
(HANE)  can  significantly  degrade  radar  target  detection  performance.  Depending  on  the 
radar  design  parameters,  similar,  but  less  severe  performance  degradation  may  occur  in 
ambient  auroral  and  equatorial  ionospheres. 
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SECTION  1 

INTRODUCTION 


A  constellation  of  space-based  radars  (SBRs)  has  been  under  consideration  eis  a  candidate 
for  the  broad  area  surveillance  and  target  tracking  required  in  future  air  defense  missions. 
These  missions  require  the  detection  and  tracking  of  aircraft  and/or  cruise  missiles  that 
could  threaten  naval  battle  groups  and  the  continental  United  States  (CONUS).  Several 
different  combinations  of  missions  amd  threats  have  been  considered  in  system  architecture 
studies.  Studies  of  the  fleet  defense  mission  consider  long  range  bombers  as  the  primary 
threat  whereas  studies  of  strategic  air  defense  of  CONUS  consider  a  combined  threat  con¬ 
sisting  of  bombers,  air-launched  cruise  missiles  (ALCMs)  and  sea-launched  cruise  missiles 
(SLCMs)., 

The  radar  propagation  path  must  traverse  the  ionosphere.  It  is  known  that  ionospheric 
disturbances  can  cause  disruption  of  transionospheric  radio  wave  propagation  (Refs.  1-16). 
Ambient  high-altitude  and  equatorial  region  ionospheric  disturbances  often  cause  severe 
amplitude  and  phase  scintillation  of  VHF  signals  (Refs.  2-5  and  7-11).  Scintillation  of 
signals  at  frequencies  up  to  L  band  is  commonly  observed  in  the  equatorial  region  due  to  a 
severe  ionospheric  disturbance  known  as  equatorial  spread  F  (Refs.,  2,4-7  and  12-14).  Phase 
scintillation  of  L-band  signals  is  also  occeisionally  observed  in  high-latitude  regions  due  to 
naturally  occurring  disturbances  (Ref.  3).  In  addition,  very  severe  scintillation  will  result 
at  frequencies  up  to  EHF  for  propagation  through  the  disturbed  ionospheric  environment 
occurring  following  a  high-altitude  nuclear  explosion  (HANE)  (Refs.  17-26). 

Ionospheric  scintillation  effects  decrease  with  decreasing  radar  wavelength,  whereas  the 
cross  section  of  cruise  missiles  tend  to  increase  with  increasing  wa,velength.  Although  a 
UHF  space-based  radar  has  been  proposed  as  a  compromise  between  these  two  trends, 
this  paper  analyzes  the  effect  of  a  disturbed  propagation  environment  on  VHF,  UHF,  and 
L-band  surveillance  radar  systems. 

Since  the  antenna  footprint  is  many  orders  of  magnitude  larger  than  the  size  of  targets,  the 
backscattered  signal  will  b  *  dominated  by  the  contribution  from  surface  clutter.  Doppler 
processing  to  discriminate  between  targets  and  clutter  is  complicated  by  the  large  Doppler 
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shift  of  the  clutter  signal  due  to  the  large  platform  velocity  of  a  SBR.  Use  of  a  displaced 
phcise  center  antenna  (DPCA)  may  be  required  at  long  radar  wavelengths  \^due  to  the  ex¬ 
tremely  large  antenna  footprints)  in  order  to  achieve  the  level  of  clutter  rejection  necessary 
for  target  detection. 

Propagation  effects  can  degrade  the  performance  of  a  space-based  surveillance  radar  in 
several  ways  (Refs.  27-30).  Amplitude  fading  can  reduce  the  probability  of  detection. 
In  addition,  phase  scintillation  can  spread  the  Doppler  spectrum  of  a  target  return  signal 
beyond  the  Doppler  resolution  of  the  radar;  this  reduces  the  ratio  of  target  signal  to  clutter 
signal  and/or  noise  with  a  resulting  reduction  in  probability  of  detection.  There  is  also  a 
reduction  in  radar  antenna  gain  and  an  increase  in  beamwidth  caused  by  a  loss  of  signal 
coherence  along  the  antenna  aperture.  Finally,  propagation  effects  can  cause  a  marked 
reduction  in  the  clutter  rejection  performance  of  the  DPCA. 

The  report  is  organized  as  follows.  Section  2  gives  a  brief  qualitative  description  of  the 
operation  of  the  DPCA,  derives  the  radar  equation  for  detection  of  a  target  in  clutter,  and 
demonstrates  why  the  use  of  a  DPCA  is  often  proposed  for  a  low-altitude  UHF  surveillance 
radau:.  The  derivation  of  the  radar  equation  includes  the  effects  of  propagation  through  the 
disturbed  environment  following  a  HANE  event.  Section  3  gives  examples  of  the  effects 
of  disturbed  propagation  environments  upon  low-altitude  UHF  and  L-band  surveillance 
radars,  while  Section  4  contains  some  concluding  remarks.  An  analysis  of  propagation  ef¬ 
fects  on  DPCA  clutter  rejection  performance  and  scintillation-induced  Doppler  broadening 
of  the  radar  signal  is  given  in  Appendix  A.  Details  of  the  propagation  environment  are 
given  in  Appendix  B. 
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SECTION  2 

INTRODUCTION  TO  THE  MODEL 


This  section  introdvices  the  model  and  illustrates  the  operation  of  pulse-Doppler  space- 
based  surveillance  radars  through  a  derivation  of  the  radar  equation  for  target  detection. 
In  order  to  make  the  discussion  of  this  section  more  concrete,  specific  parameters  for  a  UHF 
SBR  are  used  in  the  illustrations.  The  formalism  for  treating  the  impact  of  propagation 
effects  on  the  signal-to-clutter  ratio  is  introduced. 

Targets  must  be  detected  in  the  presence  of  large  surface  clutter  returns.  An  equation 
for  the  target  signal-to-clutter  power  ratio  is  derived;  this  radar  equation  includes  the 
effects  of  propagation  through  a  disturbed  ionospheric  environment.  It  is  demonstrated 
that  a  displaced  phase  center  antenna  (DPCA)  is  often  recommended  for  a  low-altitude 
UHF  surveillance  radar  under  undisturbed  propagation  conditions  to  reduce  the  clutter 
amplitude  to  acceptable  levels. 

The  effect  of  propagation  through  a  disturbed  ionospheric  environment  enters  the  radar 
equation  in  several  ways;  the  DPCA  clutter  rejection  performance  is  reduced,  the  target 
signal  is  Doppler  spread  and/or  Doppler  shifted,  the  antenna  gain  is  reduced,  and  the  an¬ 
tenna  beam  is  broadened  resulting  in  a  larger  clutter  cross  section  and  Doppler  bandwidth. 

The  displaced  phase  center  antenna  (DPCA)  is  used  to  reject  the  signal  due  to  stationary 
surface  clutter.  This  is  accomplished  by  dividing  the  antenna  into  several  subapertures 
with  phase  centers  that  are  physically  separated.  Through  a  combination  of  space  and 
time  diversity,  the  return  signal  from  surface  clutter  is  rejected.  This  technique  enables 
detection  of  moving  targets  with  radar  returns  that  are  several  orders  of  magnitude  weaher 
than  surface  clutter. 

The  simplest  implementation  of  DPCA,  which  we  call  two  phase  center  DPCA,  is  shown 
in  Figure  1.  There,  for  simplicity,  we  show  two  separate  anten’.as  whereas  actually  a  single 
antenna  is  used  with  phase  centers  that  are  electronically  adjustable.  The  operation  is  as 
follows.  At  some  time,  say  <  —  0,  a  pulse  is  tiansiuitted  from  the  front  antenna.  This 
signal  propagates  to  the  targets  and  clutter,  is  scattered,  and  propagates  back  to  the  radar 
platform.  For  a  space-based  system  the  signal  round  trip  travel  time  <</  can  be  reasonably 
large,  so  the  platform  can  move  a  considerable  distance  as  shown  in  Figure  1.  The  scattered 
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signal  is  received  by  the  front  antenna.  Due  to  platform  motion,  at  some  time  <  =  A  the 
rear  antenna  has  moved  into  the  same  position  that  was  occupied  by  the  front  antenna  at 
f  =  0.  (The  platform  is  now  shown  drawn  in  dashed  lines.)  Clearly,  this  time  A  is  given  by 
A  =  Vjd,  where  V  is  the  velocity  of  the  platform  and  d  is  the  separation  of  the  two  phase 
centers.  The  rear  antenna  now  emits  a  signal  which  propagates,  scatters,  and  propagates 
back  to  the  platform  to  be  received  by  the  rear  antenna.  As  shown  in  Figure  1,  at  time 
t  +  A  the  rear  antenna  has  moved  into  the  position  that  the  front  antenna  occupied 
when  it  received  a  signal  at  time  t  =  td.  For  a  stationary  patch  of  clutter,  it  is  evident  that 
the  two  ray  paths  are  identical,  and  hence  the  clutter  contributions  to  the  two  received 
signaJs  will  be  identical.  The  clutter  contributions  of  these  two  signals  will  cancel  when 
the  signals  are  subtracted.  On  the  other  hand,  the  signal  round  trip  time  for  a  moving 
target  will  be  different  from  the  two  pulses,  hence  there  will  be  a  residual  when  the  signals 
are  subtracted. 

There  are  a  number  of  disadvantages  to  the  two  phase  center  approach  described  in  Fig¬ 
ure  1.  In  particular,  the  pulse  repetition  frequency  (PRF)  out  of  the  two  pulse  canceller 
is  one-half  the  PRF  of  the  transmitted  waveform.  In  fact,  the  usually  recommended  im¬ 
plementation  for  both  airborne  and  space-bome  DPCA  systems  is  a  three  phase  center 
approach.  The  three  phase  center  DPCA  enjoys  the  advantage  that  it  can  be  implemented 
using  the  sum  and  difference  channels  of  a  monopulse  radar  antenna. 

The  three  phase  center  implementation  is  shown  in  Figure  2.  It  operates  as  follows.  At 
time  <  =  0  a  signal  is  transmitted  employing  the  entire  aperture  with  its  phase  center  at 
the  center  of  the  antenna.  The  signal  propagates  to  the  clutter  and  targets,  is  scattered, 
and  propagates  back  lo  the  radar  platform.  Upon  reception  at  time  t  =  td,  the  antenna 
is  configured  to  have  two  phase  centers  (forward  and  rear)  that  both  receive  the  signal. 
At  time  t  =  A,  the  platform  has  moved  such  that  the  center  phase  center  is  located  at 
the  position  that  the  front  phase  center  occupied  on  the  previous  pulse.  (The  platform 
is  now  shown  drawn  in  dashed  lines.)  We  have  the  relation  A  =  d/2V,  where  d  is  the 
separation  between  the  front  and  rear  phase  centers.  This  second  pulse  is  transmitted 
and  eventually  received  at  time  <  =  -f  A  by  both  the  forward  and  rear  phase  centers. 
Between  pulses,  the  phase  centers  have  advanced  a  distance  r//?.  The  signal  received  by 
the  forward  phase  center  for  the  first  pulse  (solid  line)  and  by  the  rear  phase  center  for  the 
second  pulse  (dashed  line)  are  shown  in  Figure  2.  The  difference  signal  between  these  two 
signals  is  taken.  This  procedure  will  null  returns  from  stationary  objects.  Cancellation  of 
the  difference  signal  from  stationary  clutter  occurs  becau.se  the  two-way  propagation  path 
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Figure  1.  Two  phase  center  DPCA  operation.  The  solid  line  is  the  first  pulse  and 
the  dashed  line  is  the  second  pulse.  Platform  motion  is  to  the  right. 

lengths  coincide  to  a  high  degree  of  accuracy.  Alternatively,  the  two-way  propagation  path 
lengths  to  a  moving  target  will  not  coincide. 

Clutter  cancellation  will  be  impaired  if  the  characteristics  of  the  DPCA  propagation  paths 
vary  between  radar  pulses  as  can  occur  if  the  radar  line-of-sight  p^lsses  through  struc¬ 
tured  ionospheric  plasma.  The  two  implementations  of  DPCA  are  affected  differently  by 
propagation  effects.  Although  the  propagation  paths  used  in  the  two  phase  center  DPCA 
physically  coincide,  they  are  separated  in  time  by  the  amount  A.  Hence,  a  time  evolution 
of  the  ionospheric  plasma,  as  occurs  when  there  is  a  bulk  plasma  drift  across  the  line-of- 
sight,  can  introduce  diff'erent  phase  shifts  into  the  two  propagation  i)aths.  Even  if  A  is  a 
small  fraction  of  the  signal  decorrelation  time,  there  can  be  a  substantial  reduction  in  the 
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Figure  2.  Three  phase  center  DPCA  operation:  The  solid  line  is  the  first  pulse 

and  the  dashed  line  is  the  second  pulse.  Platform  motion  is  to  the  right. 


clutter  rejection  performance,  with  a  consequent  reduction  in  signal-to-cluttcr  ratio.  The 
propagation  paths  used  in  the  three  phase  center  DPCA  differ  in  both  space  and  time. 
Hence,  structured  ionospheric  plasma  can  degrade  DPCA  performance  even  in  the  absence 
of  a  bulk  plasma  drift.  If  the  path  separation  is  a  small  fraction  of  the  signal  decorrelation 
length,  a  substantial  reduction  in  clutter  rejection  performance  can  result. 

Although  the  following  radar  equation  is  applicable  to  any  space-based  surveillance  radar, 
for  purposes  of  illustration,  specific  numbers  will  be  evaluated  for  a  generic  UHF  radar 
operating  at  430  MHz,  in  a  355  nautical  mile  (650  km)  circular  orbit  with  an  antenna 
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apertxire  measuring  Dax  =  48  m  in  the  azimuth  dimension  (the  antenna  dimension  par¬ 
allel  to  the  platform  motion  direction)  and  Dgi  =  16  m  in  the  elevation  dimension  (the 
antenna  dimension  perpendicular  to  the  platform  motion  direction).  The  radar  waveform 
is  assumed  to  have  a  bandwidth  of  1  MHz  and  a  pulse-repetition-frequency  of  2000  Hz. 
Calculations  will  be  performed  for  a  line-of-sight  grazing  angle  of  a  =  25°.  This  is  one  of 
two  systems  for  which  more  detailed  results  are  given  in  Section  3. 


The  resulting  antenna  footprint  measures 


AA  = 


TT  ]  ^  i?oA 

2en2\  2D  ax' 


in  the  azimuth  or  along-track  dimension  by 

^  iZoA 

2D  si  sino 


AR  = 


TT 


2£n2 


(1) 


(2) 


in  the  range  direction.  In  Eqs.  1  and  2,  Eg  is  the  slant  range,  A  is  the  radar  signal  wave¬ 
length,  and  a  is  the  radar  line-of-sight  grazing  angle.  The  footprint  widths  are  computed 
by  taking  the  3-dB  beamwidth  (in  the  absence  of  propagation  effects)  to  be  RgX/D,  fitting 
this  3-dB  beamwidth  to  a  Gaussian,  and  integrating  over  the  two-way  Gaussian  beamwidth 
to  find  the  effective  beamwidth.  The  slant  range  is  given  by 

Ro  =  [/Zg  sin^  a  +  +  2/iiZe]  ^  —  Re  sin  a  (3) 


where  h  is  the  radar  platform  altitude  and  Re  is  the  radius  of  the  Earth.  For  the  (JHF 
radar  mentioned  above,  /?<,  =  1300  km,  AA  =  14  km  and  AiZ  =  100  km  at  a  grazing  angle 
of  25°. 


As  depicted  in  Figure  3,  at  small  grazing  angles  such  that  the  range  ambiguity  length  iZamb 
is  less  than  AiZ,  surface  clutter  is  received  simultaneously  from  several  ambiguous  range 
swaths;  the  sej^aration  Ramb  of  the  ambiguoxis  range  swaths  depends  on  the  interpulse 
period  A  through 

„  e  A  . , . 

iZamb  —  X  )  (4) 

2  COS  a 

where  r  is  the  velocity  of  liglit.  The  main-beam  clutter  cross  section  <7e  equals  the  clutter 
cross  section  p^'r  unit  aiea,  ag,  muitipiied  by  the  area  of  the  ambiguous  range  swaths  lying 
within  the  antenna  footprint.  Thus  we  have 

Oc  =  naoAAArj,  (5) 
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Figure  3.  Antenna  footprint  with  ambiguous  range  swaths  hatched. 


where  n  is  the  number  of  ambiguous  range  swaths  and  the  ground  range  resolution  Ay  is 
given  by 


‘-•y  —  oE> 

2B  CCS  a 

In  Eq.  6,  B  is  the  radar  bandwidth  (assumed  to  be  1  MHz  for  the  UHF  example), 
small  grazing  angle  limit  where  Ramb  ^  ^Ri  that  is  for 


(6) 
In  the 


tano  <C 


TT  ^  iZoA 
2in2  cDeiA' 


(7) 


the  number  of  ambiguous  range  swaths  is  given  approximately  by  AR/Ramh-  Thus,  for 
Ramb  <  ^R  the  surface  clutter  cross  section  is  approximately 


TT  aoX^Rl 
8bi2  DazDefA  sin  o 


(8) 


Alternatively,  when  Ramb  >  AR,  there  can  be  at  most  one  clutter  swath  of  width  Ay  so 
that 

.[  ^  cTcXcRo 

^  2tn2\  ADazBcosa 

A  representative  value  of  <To  for  land  clutter  at  UHF,  that  will  be  exceeded  10  i)errent  of 
the  time,  is  —33  dB  (Ref.  31).  For  a  worst  case  we  will  take  (7„  =  —20  dB.  For  the  UHF 


radar  example,  Eq.  7  implies  that  the  clutter  cross  section  is  given  approximately  by  Eq.  9. 
Substitution  of  the  UHF  ladar  parameters  into  Eq.  9  and  evaluating  for  a  grazing  angle 
of  25®  gives 

<rc  =  43.7  dBm^  (10) 

However,  the  clutter  magnitude  is  significantly  reduced  through  a  combination  of  three 
methods;  Doppler  processing,  antenna  pattern  rejection,  and  use  of  a  DPCA. 

As  discussed  previously,  DPCA  will  be  used  to  cancel  the  surface  clutter  signal.  The 
amplitude  and  phase  of  the  two  signals  that  are  differenced  must  be  controlled  to  a  high 
degree  of  accuracy.  In  order  to  maintain  40  dB  of  clutter  cancellation  the  amplitude  must 
be  maintained  to  within  1  percent  and  the  standard  deviation  of  the  phase  error  must 
be  maintained  to  within  0.5  degrees.  Furthermore,  this  accuracy  must  be  maintained 
over  the  entire  antenna  pattern.  Implementation  of  these  hardware  tolerances  is  ''jry 
difficult  to  achieve,  particularly  in  the  sidelobes  of  the  antenna.  Hence,  even  in  the  absence 
of  propagation  effects,  hardware  tolerances  will  probably  limit  clutter  cancellation  to  no 
better  them  40  dB. 


The  DPCA  residual  clutter  spectrum  is  schematically  illustrated  in  Figure  4.  Although 
ideal  operation  of  the  DPCA  makes  the  radar  platform  look  stationary,  the  residual  clutter 
spectrum  due  to  imperfect  cancellation  will  have  a  spectrum  similar  to  a  non-DPCA  an¬ 
tenna.  The  most  obvious  difference  is  that  the  residual  sidelobe  spectrum  will  be  somewhat 
high  relative  to  the  residual  mainlobe  clutter  due  to  difficulty  in  matching  antenna  side- 
lobes.  The  sidelobe  clutter  occurs  between  — 2V/A  and  -{-2Vj\  in  Doppler  frequency  (where 
is  *he  platform  velocity)  while  the  mainlobe  clutter  residual  is  centered  at  {2V jX)  cos  6a 
whei'e  6a  is  the  angle  between  the  antenna  bore-sight  and  the  platform  velocity  '^ector. 
(For  the  example  UHF  radar  under  discussion,  6a  is  taken  to  be  90°. )  When  the  antenna 
is  pointed  broadside,  as  in  Figure  4,  the  mainlobe  clutter  residual  is  centered  about  zero 
frequency.  In  the  general  ca.se  there  is  also  a  large  sidelobe  return  at  zero  frequency.  This 
return,  known  as  the  altitude  return,  is  due  to  the  large  specular  clutter  return  coming 
from  the  region  directly  below  the  platform.  It  is  masked  by  the  mainlobe  clutter  when  the 


antenna  is  pointed  broadside.  The  width  of  the  inainlobe  clutter  residual  A  fd  i.s  (assuniing 


a  Gaussian  antenna  patten,  see  Appendix  Aj 


A/c/  = 


r  »  d  r 


(11) 
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The  spectrum  shown  in  Figure  4  is  not  strictly  correct  as  it  is  the  spectrum  appropriate  for 
a  CW  radar.  In  reality,  the  radar  signal  is  modulated  with  a  signal  that  typically  consists 
of  a  pulse  waveform  that  is  repeated  at  one  or  more  PRFs  over  a  finite  dwell  time  Tc-  The 
spectrum  of  a  typical  modulation  signal  is  shown  in  Figure  5.  The  spectrum  consists  of 
narrow  spikes  at  the  harmonics  of  the  PRFs.  The  spikes  have  a  width  that  is  consistent 
with  the  attainable  Doppler  resolution  (1/Tc).  The  amplitude  of  the  spikes  is  modulated 
by  the  spectrum  of  a  single  pulse;  hence,  the  overall  bandwidth  is  the  bandwidth  of  a  single 
pulse. 


DOPPLER  FREQUENCY 


Figure  4.  Doppler  spectrum  for  CW  radar. 


The  actual  Doppler  spectrum  is  the  convolution  of  the  modulation  spectrum.  Figure  5,  with 
the  spectmm  for  a  CW  radar.  Figure  4.  The  primary  effect  of  this  convolution  process  is 
that  the  spectrum  shown  in  Figure  4  is  repeated  at  the  harmonics  of  the  PRF.  With  this 
in  mind  we  return  to  our  discussion  of  Figure  4. 

Two  target  returns  are  shown  in  Figure  4  in  addition  to  the  residual  clutter  signal.  In  this 
example  one  target  is  in  the  mainlobe  clutter  region  and  another  target  is  in  the  sidelobe 
clutter  region.  Since  the  sidelobe  clutter  is  highly  Doppler  ambiguous  for  any  practical 
PRF,  sidelobe  clutter  will  occupy  all  portions  of  the  Doppler  spectrum. 

The  radar  uses  two  techniques  to  discriminate  between  various  targets  and  clutter.  At 
any  given  instant  in  time,  as  shown  in  Figure  3,  the  radar  illuminates  numerous  range 
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Figure  5.  Modulation  spectrum  of  pulse-DoppIer  radar. 

ambiguous  swaths  in  the  antenna  footprint.  The  received  signal  prior  to  range  gating  will 
contain  contributions  from  all  portions  of  the  antenna  footprint.  Range  gating  selects  one 
particular  set  of  range  ambiguous  swaths.  Hen^e,  some  targets  may  be  isolated  from  others 
through  range  gating.  The  other  discrimination  technique  is  Doppler  processing. 


As  illustrated  in  Figure  4,  the  clutter  signal  has  a  Doppler  bandwidth  Afd  t.hat  is  con¬ 
siderably  larger  than  the  Doppler  resolution  Afp  obtained  through  coherent  processing; 
thus,  only  a  small  fraction  of  the  total  clutter  signal  energy  will  lie  in  a  given  Doppler  bin. 
Coherent  processing  over  a  dwell  time  Tc  yields  a  Doppler  processing  bandwidth  Afp  of 

Afp  =  1/Tc.  (12) 


An  improvement  factor  due  to  Doppler  processing  Ip 


call  ue  ucimcu  as 


Ip  = 


A  hi 
Afp 


TT 


2Cn2 


Dn  r 


sin  6  a  ■ 


(13) 


Returning  to  the  UHF  radar  example  and  using  a  typical  coherent  processing  time  of 
=  0.1  sec  and  a  platform  velocity  of  F  =  7.5  km/s,  we  find  Afci  =  235  Hz,  Afp  --  10  Hz, 
and  Ip  =  13.7  dB. 
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Due  to  target  motion,  the  Doppler  shift  of  a  target  will  usually  differ  from  the  center 
frequency  of  the  clutter  spectrum.  As  shown  in  Figure  6,  if  the  antenna  beamwidth  is 
sufficiently  narrow  in  azimuth,  some  clutter  rejection  will  be  obtained  through  antenna 
pattern  rejection.  This  antenna  pattern  improvement  factor  I  a  can  be  estimated,  assuming 
a  Gaussian  antenna  beam  pattern  and  a  sufficiently  small  target  Doppler  shift,  through 
the  factor  (see  Appendix  A) 


I  A  =  exp 


»(2<.,/A)= 


A/. 


=  exp 


cl 


8en2v{Dl, 

V2\2 


(14) 


which  represents  the  relative  clutter  amplitude  evaluated  at  the  Doppler  frequency  corre¬ 
sponding  to  a  target  with  a  velocity  component  edong  the  radar  line-of-sight  of  Vf.  (Again, 
it  is  assumed  that  the  antenna  is  pointing  broadside.)  Assuming  a  desired  minimum- 
detectable  velocity  of  45  m/s  and  using  this  value  for  Vt  in  Eq.  14,  gives  I  a  =4.1  dB  for 
the  parameters  of  the  UHF  ratdar. 
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MAIN  LOBE  CLUTTER 


Figure  6.  Clutter  spectrum  befon*  and  after  Doppler  filtering. 


A  radar  equa^^ion  for  the  signal-to-clutter  ratio,  S/C,  can  be  written  as 


c/r^_  {^i)h^lAh)PCA 

(7cLb 


(15) 
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whore  (cr^)  is  the  average  target  cross  section,  Idpca  is  the  DPCA  improvement  factor 
and  Lb,  Eq.  A65,  is  a  target  beam  loss  factor  that  accounts  for  targets  that  are  not  in 
the  center  of  the  ant^'iina  footprint.  The  signal-to-clutter  ratio  must  be  on  the  order  of 
10  to  15  dB  for  reliable  target  detection;  the  exact  value  will  depend  on  the  probabilities 
of  detection  and  false  alarm  required,  the  statistics  of  cross  section  fluctuations,  and  the 
amount  of  incoherent  averaging. 

For  the  UHF  system,  requiring  a  signal-to-clutter  ratio  of  at  least  13  dB  implies  that  50  dB 
of  DPCA  improvement  is  required  for  detection  of  a  0.1  m^  cro.ss  section  target  (assuming 
a  minimum  detectable  velocity  of  45  m/sec).  Similarly,  40  dB  of  DPCA  improvement  is 
required  for  detection  of  a  1-m^  cro.ss  section  target.  Noting  that  the  DPCA  improvement 
factor  depends  on  target  velocity  (—5  dB  for  this  example,  Eq.  A67a),  the  clutter  rejections 
reqviired  are  55  dB  and  45  dB,  respectively. 

Thus,  good  cliiJter  rejection  performance  is  a  necessity  for  a  UHF  space-based  surveil¬ 
lance  radar  system.  The  majority  of  clutter  rejection  is  provided  by  the  DPCA.  System 
mismatches,  such  as  imperfectly  matched  antenna  patterns  from  the  two  displaced  phase 
centers,  place  a  finite  limit  on  the  cluttet  rejection  performance  of  DPCA.  Ionospheric 
induced  proj)agation  disturbances  will  further  degrade  DPCA  performance  as  well  as  in¬ 
troduce  a  scintillation-induced  loss  caused  by  the  Doppler  spread  of  target  returns.  The 
remainder  of  this  section  introduces  these  effects  into  the  radar  equation  for  the  signal-to- 
clutter  ratio. 


As  shown  in  Appendix  A,  propagation  disturbances  place  an  upper  limit  on  DPCA  perfor¬ 
mance.  The  best  achievable  DPCA  improvement  factor  Idpca  for  situation  in  which  the 
pri>pagation  path  passes  through  an  ionospheric  propagation  disturbance,  reprinted  here 
from  Appendix  A,  is  /dpca  =  Ar/Ac,  where  the  difference  signal  for  the  target  At  is 


At  =  a  —  2B  cos  2a, 

and  the  difference  signal  for  the  clutter  Ac  is 


A  i  rv 

^ 


r 

1  ^ 

rF?Ai 
1  1 

Q 
_ 1 

iO  exp  -  ^ 

For  the  two  phase  center  implementation  the  A  and  D  coefficients  are 

A2  =  2{l  +  vxp{-2ilvySa.(l)) 


(10) 


(17) 


(18a) 
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and 


n  ^  A  2A2  2 

B2  =  ::Mexp--—^V2  . 


(18b) 


The  antenna  scatter  coefficient  Sax  represents  the  reduction  in  antenna  gain  resulting  from 
propagation  effects  and  is  defined  through 


Sax  =  1  + 


4£n2  Dl, 

V  ei  ’ 


in  which  (0  is  the  decorrelation  length  for  one-way  propagation  evaluated  at  the  radar 
platform  altitude. 

For  the  three  center  implementation  the  A  and  B  coefficients  are 

Ai  =  2  +  exp[-2{tdV  -f-  d/2)‘^ / Sax^l]  +exp[-2(tdV  -  d/2)‘^ / Sax^l]  (20a) 


where 


i?3  =  exp  ~  [exp{-d'^ /2Sax(l)  +  exp(-2t^F^/5a-^o)]  (20b) 


d  =  2A{V-VEx).  (20c) 


The  following  symbols  have  also  been  defined: 

d  =  —ARq  •  Vo  -1 - ^ +  -~VEx  +  , 

V2  =  VEx  +  (21) 

We  have  specialized  to  a  one-dimensional  phase  screen  with  variation  and  a  possible  drift 
motion  along  the  platform  motion  direction  (x  direction).  is  the  distance  between  the 
platform  and  the  pha.se  screen  and  i?-  is  the  distance  between  the  phase  screen  and  the 
earth’s  surface.  Rq  is  a  unit  vector  pointing  from  the  antenna  center  to  the  target,  V 
is  the  platform  velocity  relative  to  the  phase  screen  (which  is  by  assumption  constrained 
to  be  in  the  x  direction),  vr  is  the  earth’s  velocity  relative  to  the  pha.se  screen.  Vn  is 
the  target  velocity  relative  to  the  earth,  and  Vx  represents  either  the  x  component  of  the 
earth  velocity  (when  clutter  is  being  considered)  measured  relative  to  the  phase  screen, 
or  the  X  component  of  the  target  velocity  measured  relative  to  the  phase  screen.  Hence, 
Vx  represents  vex  when  clutter  is  being  considered  and  vex  +  when  the  target  is 
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being  considered.  Furthermore,  we  have  specialized  to  the  case  where  the  antenna  points 
broadside. 


In  addition  to  the  degradation  of  DPCA  performance,  propagation  induced  scintillation 
causes  an  apparent  Doppler  spread  and  Doppler  shift  of  the  return  signal.  Although  this 
effect  occurs  to  both  the  surface  clutter  return  and  the  return  from  a  discrete  target, 
the  propagation  induced  Doppler  spread  and  shift  of  the  clutter  is  generally  negligibly 
small  compared  to  the  main-beam  clutter  bandwidth  Afd.  Alternatively,  the  Doppler 
shift  and/or  spread  can  be  significant  compared  to  the  radar  Doppler  resolution  A/d;  the 
Doppler  spread  can  spread  the  energy  returned  from  a  discrete  target  over  a  spectral  width 
greater  than  A fo  with  a  consequent  decrease  in  signal-to-clutter  ratio.  In  addition,  a  net 
propagation  induced  Doppler  shift  may  cause  the  target  to  appear  in  the  wrong  Doppler 
bin.  This  maj"^  not  be  important  if  accurate  target  velocity  information  is  not  needed,  as 
may  be  the  case  when  the  radar  is  in  the  search  mode.  Alternatively,  if  the  net  Doppler 
shift  changes  from  scan  to  scan,  and  if  some  incoherent  processing  is  performed,  there  may 
be  a  further  performance  reduction. 


The  performance  reduction  due  to  a  Doppler  spread  is  quantified  through  a  coherent 
integration  loss  factor  Tc«-  When  the  scintillation  induced  Doppler  spread  Af$c  is  less 
than  the  radar  Doppler  resolution  Afo,  there  is  no  loss. 


Lc,  =  1. 


(22) 


When  A/,c  >  A//),  the  coherent  integration  loss  factor  becomes 

Lc.=  ^  =  AUT,.  (23) 

If  the  Doppler  spread  is  due  primarily  to  multipath  propagation,  that  is,  when  the  contri¬ 
bution  due  to  the  phase-only  effects  can  be  neglected,  the  coherent  integration  loss  trom 
Appendix  A  becomes 


L 


CA  — 


,  2 

\2V^T^] 

1  +  - 

TT 

2V^T^ 
H - - — ^ 


1 

2 


(24) 


We  assume  that  the  radar  system  will  be  designed  such  that  under  normal  operating 
conditions,  that  is,  in  the  absence  of  ifjno.spheric  induced  propagation  distuibances,  the 
noise  level  will  be  greater  than  the  clutter  signal  level  at  all  Doppler  frequencies  of  interest. 
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When  a  propagation  disturbance  occurs,  the  clutter  signal  level  will  increase  to  levels 
exceeding  the  noise  level;  under  these  disturbed  conditions  the  detection  performance  is 
given  in  terms  of  the  radar  equation  for  the  signal-to-clutter  ratio  S/C  (Eq.  A60), 

5  _  (<T()/d^/i/dpca 

C  “  <TcI=i/.i,v/K7  ■  '  ’ 

Propagation  effects  will  be  investigated  for  candidate  UHF  ?ind  L-band  systems  in  Sec¬ 
tion  3. 
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SECTION  3 

PROPAGATION  EFFECTS  ON  CANDIDATE  SBR  SYSTEMS 


This  section  applies  the  space-based  radar  target  detection  model  to  candidate  UHF  and 
L-band  pulse-Doppler  surveillance  radars.  The  degradation  in  signal-to-clutter  ratio  is 
quantified  for  propagation  envirorunents  ranging  from  benign  natural  to  post  HANE  envi¬ 
ronments.  Performance  of  two  phase  center  and  three  phase  center  DPCA  implementations 
are  contrasted. 

The  parameters  of  the  ceindidate  L-band  and  UHF  systems  are  listed  in  Table  1.  These 
candidate  systems  were  analyzed  primarily  to  illustrate  the  manner  in  which  the  model 
is  utilized  and  uncover  the  relative  sensitivity  to  propagation  effects;  hence  the  choice  of 
SBR  parameters  may  not  represent  the  .systems  currently  under  consideration  for  future 
deployment. 


Table  1.  Parameters  of  the  candidate  space-based  radar  systems.  The  antenna 
aperture  dimensions  are  listed  with  the  along-track  (or  azimuth) 
dimension  proceeding  the  cross-track  (or  elevation)  dimension. 


Parameter 

L-band  System 

UHF  System 

Frequency 

1.2  GHz 

430  MHz 

Altitude 

650  km 

650  km 

PRF 

5000  Hz 

2000  Hz 

Grazing  Angle 

25° 

25° 

Antenna  Aperture 

40  m  X  15  m 

48  m  X  16  m 

Radar  Bandwidth 

1  MHz 

1  MHz 

Integration  Time 

0.1  sec 

0.1  sec 

Target  Cross  Section 

0  dBsm 

0  dBsm 

Clutter  Reflectivity 

-20  dB 

-20  dB 
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An  ionospheric  striation  drift  velocity  of  200  m/sec  was  assumed.  Although  the  (circular) 
orbital  velocity  is  7.5  km/sec  for  the  selected  altitude,  all  velocities  used  in  Section  2  and 
Appendix  A  eure  measured  relative  to  the  rest  frame  of  the  striations;  hence,  a  platform 
velocity  of  V  =  7.3  km/sec  was  taken.  Depending  on  the  orbital  parameters  and  the  radar 
look  geometry,  the  velocity  of  the  Earth’s  stirface  can  range  up  to  several  hundred  meters 
per  second.  A  typical  value  (relative  to  the  ionospheric  striations)  of  vgj;  =  200  m/sec  \vas 
taken. 

Signal-to-clutter  ratios  depend  strongly  upon  the  target’s  velocity  component  along  the 
radar  line-of-sight.  All  calculations  shown  in  this  section  were  performed  for  a  typical  min 
imum  detectable  target  velocity  of  vt  =  45  m/sec.  Signal-to-clutter  ratios  were  calculated 
versus  decorrelation  length  evaluated  at  the  radar  platform  altitude.  The  phase  screen 
was  taken  to  be  located  halfway  between  the  Earth’s  surface  and  the  platform  altitude. 
(Details  of  the  propagation  model  are  given  in  Appendix  B.) 

Figure  7  shows  the  signal-to-clutter  ratio  for  the  L-band  system.  The  upper  curve  is  for 
the  two  phase  center  implementation  and  the  lower  curve  is  for  the  three  phase  center 
implementation.  As  discussed  in  Section  2,  DPCA  performance,  and  consequently  the 
signal-to-clutter  ratio,  is  limited  by  niunerous  effects  in  addition  to  propagation  effects. 
Since  these  effects  are  not  explicitly  included  in  the  model,  the  DPCA  improvement  factor 
(shown  in  Figure  9)  increases  without  bound  as  the  decorrelation  length  increases.  In 
reality,  the  DPCA  improvement  factor  saturates  at  a  finite  value,  typically  on  the  order  of 
40  dB. 

Figures  8,  9  and  10  demonstrate  the  behavior  of  the  key  elements  that  enter  the  signal-to- 
clutter  calculation.  As  these  results  show,  the  performance  degradation  due  to  propagation 
effects  results  mainly  from  a  combination  of  a  loss  of  Doppler  processing  gain  (Figiu'e  8)  and 
a  degradation  in  DPCA  clutter  rejection  (Figure  9 ).  The  DPCA  target  response  (Figure  10) 
shows  little  degradation.  The  loss  of  Doppler  processing  gain  is  due  to  the  combined  effects 
of  amplitude  fading  and  phase  scintillation.  Amplitude  fading  can  reduce  the  energy  in  the 
target  signal  and  phase  scintillation  can  spread  the  Doppler  spectrum  of  a  target  return 
signal  beyond  the  Doppler  resolution  of  the  radar;  both  show  up  as  a  reduction  in  Doppler 
processing  gain.  Losses  due  to  DPCA  clutter  rejection  occur  from  phase  scintillation  over 
the  time  separation  A  or  spatial  separation  d  of  the  two  pro})agation  paths.  A  phase 
imbalance  of  a  fraction  of  a  radian  is  sufficient  to  disrupt  clutter  rejection;  hence,  DPCA 
performance  is  disrupted  when  A  is  a  small  fraction  of  the  signal  decorrelation  time  or 
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(for  the  three  phzise  center  implementation)  when  c?  is  a  small  fraction  of  the  decorrelation 
length  Cq.  Although  severe  scintillation  can  also  cause  a  reduction  in  the  signal-to-clutter 
ratio  through  antenna  scattering  loss,  this  effect  becomes  significant  only  at  extremely  small 
decorrelation  lengths  (at  which  Doppler  processing  losses  and  DPCA  losses  are  severe). 


Figure  7.  Signal-to-Clutter  ratio  as  a  function  of  decorrelation  length  for  the 
peuameters  of  Table  1:  L-band. 


As  expected,  the  two  phase  center  DPCA  implementation  is  considerably  more  robust 
against  propagation  effects  than  the  three  phase  center  implementation.  In  fact,  the  two 
phase  center  implementation  maintains  a  15  dB  signal-to-clutter  ratio  down  to  decor¬ 
relation  lengths  of  10  m,  whereas  the  signa-l-to-clutter  ratio  for  the  three  pha.se  center 
implementation  drop.s  below  15  dB  for  decorrelation  lengths  less  than  100  in. 


A  corresponding  series  of  plots  is  shown  for  the  UHF  system  in  Figures  11  14.  The 
same  general  trends  are  evident.  Examinafion  of  these  figures  reveals  that  the  candidate 
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UHF  system  could  experience  propagation-related  performance  disruptions  due  to  natu¬ 
rally  occurring  ionospheric  disturbances.  In  particular,  for  the  parameters  discussed  in 
Appendix  B  for  severe  spread  F{£o  ~  250  m),  the  two  phase  center  system  would  have 
a  signal-to-clutter  ratio  of  18  dB  while  a  value  of  0  dB  would  apply  for  the  three  phase 
center  approach. 


Figure  8.  Doppler  processing  gain  as  a  function  of  decorrelation  length:  L-band. 
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Figure  1 


SECTION  4 

CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  WORK 


A  comprehensive  model  for  the  impact  of  ionospheric  propagation  effects  on  pulse-Doppler 
space-based  rada  3  in  a  target  search  mode  has  been  developed.  Propagation-induced 
DPCA  degradation,  Doppler  processing  losses,  and  antenna  scattering  losses  are  all  prop¬ 
erly  incorporated.  Our  model  differs  from  previous  models  of  space-based  radars  with 
DPCA;  previous  models  (Refs.  27-30)  assumed  that  propagation  effects  had  a  negligible 
effect  on  DPCA  cancellation.  These  previous  studies  quantified  processing  losses  but  did 
not  quantify  DPCA  losses  due  to  propagation  effects.  We  find  that  DPCA  processing 
losses  can  be  significantly  larger  than  Doppler  processing  losses. 

Although  the  basic  model  is  applicable  to  any  degree  of  propagation  disturbance,  analytic 
results  are  possible  only  for  Rayleigh  scintillations.  Phase  screen  simtjlations  are  necessary 
for  less  severe  propagation  environments  for  whieh  non- Rayleigh  statistics  characterize  the 
one-way  path  (Ref.  32). 

Analytic  results  were  presented  for  candidate  L-bzind  and  UHF  space-based  radar  under 
Rayleigh  scintillating  conditions  ranging  from  relatively  benign  naturally  occurring  con¬ 
ditions  to  post-HANE  conditions.  All  candidate  systems  were  susceptible  to  considerable 
performance  degradation  in  the  post-HANE  environment,  although  the  two  phase  center 
li-band  system  was  surprisingly  robust.  In  all  cases  the  two  phase  center  implementation 
was  considerably  more  robust  than  the  three  phase  center  implementation.  The  UHF 
system  appears  to  be  susceptible  to  performance  degradation  resulting  from  naturally  oc¬ 
curring  ionospheric  conditions. 

There  are  several  obvious  mitigants  for  propagation  effects.  The  most  obvious  mitigant  is 
the  use  of  higher  operating  frequencies.  Another  mitigant  is  the  use  of  a  larger  antenna 
aperture  in  the  along  track  direction.  This  narrows  the  Doppler  spectrum  of  the  clutter 
and  produces  less  reliance  on  the  very  critical  DPCA  phase  cancellation. 

A  Hiwei  platform  altitude  is  another  possible  mitigant.  This  produces  less  d’fference  in  tlie 
DPCA  path  lengths  and  hence  less  disruption  in  the  DPCA  clutter  cancellation.  It  may 
also  be  possible  to  avoid  propagation  through  disturbed  regions  through  use  of  sufficiently 
low  platform  altitudes. 
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The  two  phase  center  implementation  is  considerably  more  robust  than  the  three  phcise 
center  implementation.  This  is  because  although  the  two  propagation  paths  are  separated 
in  time,  they  are  over  identical  spatial  paths.  This  implies  that  the  two  phase  center 
implementation  is  susceptible  only  to  the  bulk  striation  drift. 

There  are  numerous  possibilities  for  additional  work  in  this  area.  More  detailed  param¬ 
eter  studies,  for  example,  the  variation  of  signal-to-clutter  ratio  with  target  velocity  for 
fixed  values  of  decorrelation  length,  etc.,  can  be  performed.  Phase  screen  simulations 
could  be  performed  to  more  accurately  characterize  the  onset  of  performance  degradation. 
(As  discussed  above,  the  calculations  performed  in  this  manuscript  were  done  under  the 
assumption  of  Rayleigh  scintillations.)  In  addition,  it  may  be  desirable  to  model  addi¬ 
tional  radcir  modes  such  as  a  target  tracking  mode.  (Calculations  in  this  manuscript  were 
performed  for  the  degreidation  in  signal-to-clutter  ratio  while  in  a  target  search  mode.) 
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APPENDIX  A 

RADAR  EQUATION  DERIVATION 


We  will  first  give  a  few  more  details  of  the  radar  that  enters  the  calculation  of  the  radar 
equations.  Each  phase  center  of  the  radar  has  a  (1-D)  Gaussian  aperture  illumination 
function,  a(x),  associated  vdth  it,  appropriate  to  a  side  looking  radar.  Then,  a(x)  is 

a(x)  =  ~exp[-A:^jX^]  (Al) 

s/t: 


where  kox  is  related  to  the  aperture  size  [see  Eq.  (A5)].  The  one-way  antenna  radiation 
pattern  as  a  function  of  c,  the  angle  from  broadside,  is  then 


(A2) 


where  ko{—  27r/A)  is  th 
approximated  as 


adar  wave  number.  For  small  angles,  [ej  -C  1,  G(€)  can  be 


G(€)  =  exp 


Mf!! 

2^01. 


(A3) 


If  we  require  the  3  dB  point  to  satisfy 


G 


(A4) 


where  D  is  the  azimuthal  antenna  length  (called  D^z  in  Section  2),  then  we  have 


A-2  =  — 

2^n2D2  ■ 


(AS) 


The  radar  emits  a  coded  pulse,  which,  upon  matched  filtering,  produces  a  ground  range 
/♦’solution  of  8y.  We  will  approximate  the  point  spread  function  (PSF)  as 


PSF{y-y,)  = 


1, 

0, 


\y  -  ?/oi  < 
otherwise, 


(AG) 


Depending  on  the  operating  parameters  of  the  radar,  the  range  may  or  may  not  be  am¬ 
biguous  within  the  main  lobe  of  the  radar.  If  the  range  is  ambiguous,  contributions  from 
all  the  ambiguous  range  cells  must  be  included.  The  expressions  for  the  signals  that  wc 
will  employ  will  always  be  after  match  filtering. 
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The  background  consists  of  a  collection  of  scatterers,  located  at  positions 

We  will  assume  the  background  is  completely  incoherent  so  that  a  simple  model  has  <f)j  a 
random  variable  uniformly  distributed  between  0  and  27r  and  a  constant,  given  by 

-4j=cr^/nxny  (A7) 

where  rix{rty)  is  the  number  of  scatterers  per  unit  length  in  the  x{y)  direction. 

In  the  following,  we  assume  that  the  fading  statistics  are  Rayleigh  and  separate  the  transfer 
function  as  follows  [cf.  Eq.  (B6)] 

h  =  hsexp[i<l)i  +  i4>2]  (A8) 

where  hg  represents  the  Rayleigh  scintillation  contribution  and  +^2  due  to  phase-only 
effects.  As  discussed  in  Appendix  B,  we  divide  the  spectrum  at  a  wave  number  Kfi  such 
that  larger  than  values  contribute  to  the  Rayleigh  scintillation  and  smaller  than  Kfi 
values  contribute  to  the  phase  only  effects.  If  we  define  the  length,  L/?,  as 

Lr  =  27r/X/?,  (A9) 

the  phase-only  effects  are  here  further  divided  into  two  parts.  The  rationale  for  this  division 
is  that  there  is  an  additional  scale  of  importance,  i5,  the  distance  the  radar  travels  during 
a  dwell  time.  If  Lr  <  Ls^  we  identify  the  contributions  as 

4>\  :  Lr  <  scale  size  <  L5, 

(f)2  :  Ls  <  scale  size. 

The  (l>i  contril)ution  will  produce  a  Doppler  spread  since  it  leads  to  pliase  errors  during  the 
coherent  integration  time  that  are  approximately  averaged  while  (j>2  will  lead  to  a  Dopj>ler 
shift  since  these  ]>hases,  associated  with  long  w^avelengths,  are  essentially  constant  during 
the  dwell  time.  If  Lr  <  Lr^  we  only  have  (f>2  type  contrib\itions.  This  division  of  the 
spectrum  is  rather  arbitrary  but  since  the  dominate  effects  are  produced  by  the  Rayleigh 
scintillations,  the  distinctions  are  not  critical.  In  the  following  we  assume  the  phase- 
only  effects  only  affect  the  Doppler  processing  and  are  essentially  ccnistant  for  tlie  DPCA 
operation  [sec  Eq.  (A2G)]. 

The  displaced  phase  center  antenna  is  used  to  reject  the  signal  diu'  to  stationary  surface 
clutter.  This  is  accomplished  by  making  the  phase  center  of  the  antenna  electronically 
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adjustable  as  discussed  in  Section  2.  Here,  we  will  generalize  the  process,  instead  of 
restricting  it  to  the  two  implementations  discussed  in  Section  2.  The  operation  consists  of 
defining  two  series  of  pulses.  The  first  series  consists  of  pulses  transmitted  with  the  phase 
center  at  some  location,  designated  Xi,  as  measured  relative  to  the  center  of  the  physical 
antenna,  and  then  received  at  phase  center  location  Fi.  A  second  series  is  defined  with 
the  phase  center  located  at  X2  for  transmission  and  Y2  for  reception.  These  two  series  are 
essentially  subtracted  for  surface  clutter  rejection.  For  this  to  happen,  the  location  of  the 
phase  centers  and  the  time  separation  of  pulses  have  to  be  related,  as  discussed  below. 

The  two  phase  center  concept  is  defined  by  assuming 

Xi=Yi=  d/2,  (AlOa) 

X2  =  F2  =  -d/2.  (AlOb) 

It  is  to  be  noted  that  for  a  system  deployed  in  space,  the  clutter,  in  general,  will  be  moving 
due  to  the  earth’s  rotation  so  the  ray  paths  are  not  identical  even  for  the  clutter.  As  we 
will  see,  a  slight  altering  of  the  processing  that  was  discussrd  in  Section  2  can  reject  the 
moving  clutter  and  detect  targets  moving  relative  to  the  clutter.  The  three  phase  center 
concept  is  defined  assuming 


II 

to 

II 

0 

(Alla) 

Fi  =  d/2, 

(Allb) 

il 

1 

to 

( A.llc) 

It  is  clear  for  this  implementation  that  the  ray  paths  are  not  identiced  even  for  truly 
stationary  clutter. 

The  starting  point  is  the  exi)ressions  for  the  returned  signals  for  the  two  series.  The 
coordinate  system  employed  for  the  calc\dation  is  the  rest  frame  of  the  disturbing  medium. 
In  addition  we  assume  there  is  structure  only  in  the  along  track  direction  (x,i  =  unit 
vector).  If  the  center  of  the  antenna  is  at  .r  =  0,  at  /  =  0,  a  point  on  the  antt'iina  (given 
by  ^  or  7/  relative  to  the  phase  center)  for  the  fiist  s('ries  of  pulses  is  given  l)y 

r,  -  tiTV  4-  A'l  -I-  6  , 

7/ 1  —  ( uT  -r  f  d )F  +  j  ]  r  7/1  ■ 
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(A12a) 
(A 12b) 


Here,  we  have  defined 


n  =  pulse  number, 

T  =  interpulse  period  for  series  1  (or  2), 
r  =  position  on  transmit, 
u  =  position  on  receive, 
td  =  2Ro/c  =  time  delay, 

Ro  =  antenna  —  scatterer  slant  range, 

V  =  platform  velocity. 


For  the  second  series,  the  corresponding  positions  are 

r2  =  (nT  +  A)F  +  X2+6, 

U2  =  (nT +  td  + A)V +  Y2 +rii,  (A13) 

where  A  is  the  time  separation  between  the  two  series.  Note  therefore  that  A  is  the 
time  interval  used  for  the  DPCA  processing  and  T  is  the  time  interval  used  for  Doppler 
processing.  The  scatterers  positions  for  the  two  series  are  given  by 


wi  =  Wo  +  nTv, 
W2  =  Wi  4-  Av, 


where  Wo  is  the  initial  position  for  the  first  pulse, 


Wo  =  rr  +  -tdv 


(A14a) 

(A14b) 


(A15) 


and  v  is  the  velocity  of  the  scatterer.  In  particular,  if  the  scatterer  is  clutter,  v  =  v/p,  the 
earth’s  velocity  in  the  rest  frame  of  the  media,  and  if  it  is  a  target,  v  =  vg  +  Vo,  where 
Vo  is  the  velocity  of  the  target  relative  to  the  earth.  Under  the  usual  conditions  of  a  small 
radar  footprint  and  a  reasonably  small  dwell  time,  all  slant  ranges  c  be  approximated 
as 


(A16a) 
(AlGb) 
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R  =  i^o  +  AR, 

R  =  Ro  Ro  •  AR. 


The  returned  signals  for  the  two  series  after  matched  filtering  are 


^i(«)  =  drii  a{^i)a{r}i)h{ri,wii:)h{ui,wu) 

j 

27r 

X  expt  — [2i2o  +  Rq  •  {2nTv  —  rji  —  (Al7a) 

A 

«2(n)  =  /  d^i  diji  a{(i)a{rii)h{r2,W2i)h(u2,W2t) 

j  '' 

Onr  ^ 

X  exp  i—  [2i2o  +  Rq  ■  (2nrv  +  2Av  -  rji  -  «2i)]  •  (Al7b) 

Here,  we  have  included  point  to  point  transfer  functions  from  points  on  the  antenna  to 
points  in  the  clutter  or  target.  The  integration  over  tremsmit  and  receive  apertures  will 
lead  to  finite  aperture  effects. 

The  DPCA  processed  signal  is  the  appropriately  weighted  difference  of  sj  and  S2, 

T{n)  =  si{n)-e-'^S2{n).  (A18) 

Noting  that 

r2  +  U2  =  r,  +  Ml  +  2AV  -  (JA,  +  T,  -  X2  -  ^2),  (A19) 

we  have 

T{n)  =  'Y^Aje'‘^>  j  d(idrji  a(^i )a(r;i)  exp  [2i?o  +  .Ro  •  (SnTv  -  rii  -  uii)] 

X  [h{ri,wij:)h(ui,wii)  -  e''^h(r2,W2i)h{u2,W2T)]  (A20) 

where 

7  =  y  i?o  •  [-2A(Fi  -  v)  +  (X,  +  Fi  -  X2  -  r2)i]  -  0.  (A21) 

Clearly,  for  the  ambient  environment,  we  must  require  7  to  vanish  for  all  scattcrers  that 
constitute  the  clutter.  For  a  narrow  beam  with  zero  squint  angle  and  look  angle  d,Ro  for 
fixed  range  is  approximately 

Ro  =  N  +  (l>i  (A22) 

where  N  is  the  liemn  <lirection, 

N  =  (O..sin^,-cos^)  (A23) 
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and  (f>  is  the  angle  off  be?im  center.  The  cancellation  of  the  coefficient  of  <i>  determines  A, 
the  time  separation  of  the  two  series, 


lXi+Yi~X2-Y^ 

2  V-VEr 


(A24) 


while  the  part  independent  of  <f)  determines  /?, 


47r 

/?  =  -^Asm9vEy  (A25) 

With  the  proper  choices  of  A  and  the  DPCA  processed  signal  can  be  written  as 
r(n)  =  Y,  A,e‘'^^«e‘^e‘^expi^[2i?o  +  ^  •  (2nT(v£;  +  Vq  -  Vi)  -  i{UV  +  A,  +  Fi))] 


X  J  d^idrji  a(6)«(^i)exp-t— ilo  •  i(6 +»7i) 


X  e"'°‘hs{ri,wu)hsiuuwii)  -  e'°‘h3{r2,W2x)hs{u2,W2x)  (A26) 

where,  again,  Vq  is  the  velocity  of  the  target  relative  to  the  earth.  The  phcise  only  effects 
are  contained  in  S, 


8  =  24>2{nTV,Xj)  +  <t>x{nTV,Xj)  +  4>,  ((nT  +  td)V,Xj), 


(A27)) 


which  refers  to  both  the  antenna  position,  for  example,  nTV,  and  the  location  of  the  scat- 
terer,  Xj,  and  makes  use  of  the  fact  that  these  effects  are  associated  with  long  wavelengths. 
We  have  also  defined  a  as 


^ARo  -  Vo, 


clutter, 

target. 


(A28) 


The  DPCA  signal  is  Doppler  processed  to  detect  moving  targets.  For  the  purposes  of 
analysis,  we  will  employ  a  .simple,  unweighted,  Doppler  filter,  given  by 


I'^oi  < 


(A29a) 


(A29b) 


where  N  is  the  total  number  of  pulses  in  the  coherent  integration  time.  In  particular,  the 
quantity  we  will  investigate  is  the  average  value  of  the  square. 


<|^(V^D)r)  =  ^  E  <  r(n)r(n')  >  . 


(A30) 


n,n'=0 
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We  note  here  that  since  the  background  is  modeled  as  a  set  of  incoherent  scatterers,  we 
have  the  effective  replacement 


(A31) 


J  ] 


where  Xj  is  a  sum  over  an  azimuthal  line  of  scatterers.  As  for  the  target,  we  have 


— >  <7 - 


(A32) 


For  the  phase  only  effects,  we  have 


<  e 


>=  e^‘  [MnTV,Zj)-4>i(n'TV,zj )] 


(A33) 


where 


r  =  2(71(0)  +  2Ci(tdV)  -  2Ci  [(n  -  n')TV] 

-  Cl  [(n  -  n')TV  +  tdV]  -  Cl  [(u  -  n')TV  -  UV] .  (A34) 


Here,  C7i(x)  is  the  phase  correlation  function  that  arises  from  that  part  of  the  spectrum 
that  contribvites  to  (f)i. 

The  Rayleigh  scintillation  contribution  can  be  readily  calculated  using  Eqs.  (B9)  and  (BIO). 
In  addition,  with  the  assumed  aperture  illumination  function,  Eq.  (Al),  the  aperture 
integrals  are  straightforward.  The  fact  that  4>2  depends  on  the  scatterer  location  can 
complicate  the  sum  over  the  background  scatterers.  However,  it  is  known  that  Doppler 
processing  essentially  selects  a  restricted  region  of  the  footprint  such  that  the  scatterers 
have  the  appropriate  Doppler  shift.  Since  <^2  is  assumed  to  vary  slowly,  we  can  treat  x,  in 
Eep  (A33)  as  effectively  constant.  However,  since  the  target  and  the  clutter  it  is  competing 
against  in  a  given  Doppler  cell  are  physically  located  at  different  azimuth  positions,  the 
Doppler  shifts  arising  from  </>2  will  be  different.  (Typical  magnitudes  of  these  shifts  will 
be  discussed  in  Appendix  B.)  Therefore,  for  the  clutter  contribution,  the  s»un  over  the 


35 


background  scatterers  can  be  effectively  approximated  by  a  Gaussian  integral.  We  then 
find 


<T(n)T*(n')) 


X  exp  -  •  W) 

X  |A-2Bexp-fc^^5^^^^ 

47r^  1 


t\2 


T\n  -  n!) 


.  4A 

A2v;2cosh|^TF2^(n 

^0 


1 


Target  ■'^2^A2  5 


X  e-.-^TRo  (W-vi)(n-n')  _^y2rj.2^^  _ 

OCq 


X 


f  4A 

|a  -  Be2‘«exp-^Fil/2n«  '  «') 

4A  1 

-  Be-2‘“exp  —ViV^Tin  -  n')| 


where  (identified  as  Saz  in  Section  2) 


^  ,  2  ,  Atn2 

■5^  =  1  +  /.O.O  =  1  + 


^O'^'Or 


We  note  the  bounds 


1  1 

0  <  ^  < 


-  Sel  Un2D‘^' 


0  < 


-  ki,e„s 


<  1. 


The  effective  velocity  of  the  clutter  in  the  center  of  the  beam  is  W, 


W  =  ri-V£;-^i 

while  the  effective  target  velocity  relative  to  W  is 


T/ 


V„  =  Vo  T  — 


"')]}. 


(A35) 


(A36) 

(A37) 

(A38a) 

(A38b) 

(A39) 


r  A  AO\ 


For  DPCA,  the  target  velocity  hcis  been  effectively  changed  such  that  [cf.  Eq.  (A28)] 


n  =  — Ai?o  •  Vo  + 


5- 


-  (  VEr 


/?+  ,  , 


(A41) 
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We  have  also  defined  various  velocities  as 


Vi  =  F  +  |it;,, 

rr 

V2  =  VEx  + 


(A42) 


(A43) 


The  clutter  contribution  is  proportional  to  the  total  clutter  cross  section  for  a  range  cell, 


(Tc  =  —=a^6yRQ\ko^. 
2\/ir 


The  coefficients  A  and  B  are  given  as 


(A44) 


A  =  2  +  exp  -2(A:i  -  T,  -  tdVf/Sll  +  exp  -2(X2  -  F2  -  UYffSe^,  (A45) 

9 a2  1  f  1  r  1 

B  =  exp  "se  F'  - ^2)'  +  (>'>  2<^‘  +yi-X2- 

+  exp  (V,  -X2f+  (X,  -Y^f-  i(X,  +  K,  -  X,  -  Yif 


+  20^5'!  +  n  -  X,  -  JIC2)  +  2(UV)‘ 


(A46) 


The  particular  expressions  for  the  two  and  three  phase  center  concepts  discussed  in  Sec¬ 
tion  2  are,  for  the  two  phase  center  concept. 


A2  =  2(1-1- exp -2(tdF)V5f^), 

1  2A^ 

i?2  =  -A2exp-^F2^ 


(A47a) 

(A47b) 


and,  for  the  three  phase  center  concept, 


A3  =  2  +  exp -2^dV  +  ^  /S^l+exp-2^dV /SCl,  (A48a) 


B3  =  exp-^^V2{exp-dF |2S^l  +  exp -2{tdVf  / S^l). 

Note  that  we  have  for  the  two  implementations  (cf.  Eq.  (A24)] 

two:  A  =  d/{V-VEz), 
three:  A  =  ^d/ {V  -  vex), 

while  convenient,  but  not  necessary,  choices  for  T  are 


(A48b) 


(\49) 


two:  T  =  2  A, 
I  ee:  T  =  A. 


(A50) 
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For  the  Doppler  processing,  we  will  approximate  the  phase  only  effects  as  [cf.  Eq.  (A33)] 

hiriTV.Xj)  -  <f>2{nTV,Xj)  =  C2(xj)(n  -  n'),  (A51a) 

^  =  Cl(n-n')^  (A51b) 

(see  Appendix  B).  We  also  notice  that  the  (n  —  n')  dependence  inside  the  {  }  in  Eq.  ( A35) 
and  Eq.  (A36)  have  small  coefficients  and,  to  first  approximation,  can  be  ignored.  This 
allows  the  ionospheric  effects  on  DPCA  and  Doppler  to  be  separated.  The  Doppler  pro¬ 
cessing  for  both  the  target  and  clutter  then  have  the  general  form 

n,n'=0 

=  E  (A52) 

where 

Kl  <  5.  (A53) 

The  expressions  for  ^  and  Z  for  the  clutter  and  target  can  be  determined  from  Eqs.  (A35) 
and  (A36),  respectively,  and  Eq.  (A30). 


The  maximum  response  of  the  Doppler  processed  target  signal  is  characterized  by 


F(^  =  0,  Zt)  ~ 


1-t-  -iV 

TT 


Zt  = 


(A54) 

(A55) 


The  requirement  of  ^  0  determines  a  particular  Doppler  velocity,  T  £>,  the  Doppler  bin 

with  peak  response,  that  satisfies 


=  i?„ .  (W  -  vi)  -  ^  A  +  (A66) 

where  (  is  an  integer  (the  Doppler  velocity  may  be  ambiguous)  such  that  Eq.  (A29b)  is 
satisfied  and  Ci{T'}  is  the  appropriate  value  for  the  target  location  in  the  footprint. 


The  clutter  processed  signed  at  Doppler  velocity  Vo  is  clunacterized  by 


(A57) 
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where 


Zc  =  kl,S{i  •  WfT^  +  +  Cl 


(A58a) 

(A58b) 


(A59) 


where  t  is  an  integer  such  that  Eq.  (A53)  is  satisfied  and  C2(C)  is  the  appropriate  value 


for  the  interfering  clutter. 


The  above  simplifications,  including  the  separation  of  DPCA  and  Doppler  effects,  lead  to 
the  following  simple  approximate  expression  for  the  signal  to  clutter  ratio. 


S/C  —  ^  ^upcaIqIa 
^  “<^0  S^LBLci 


In  the  above,  we  have 


(A60) 


a  ~  average  target  cross  section 
Of.  —  average  surface  clutter  cross  section 
If  7r  1  ^  0  XRo 

=  2  [itejJ 

n  =  number  of  ambiguous  range  cells 
^DPCA  =  DPCA  improvement  factor 
=  At/ Ac 

Id  Doppler  improvement  factor  for  clutter 

_  ■  F  iVTc 
[2/>n2\  D 

I A  =  antenna  pattern  improvement  factor 


=  exp  Sin2 


A2y2 


Lfl  =  beam  loss  due  to  target  location(72o) 

D'^  1  - 

=  exp 8£n2—-(Ro  •  i) 

Lfi  =  coherent  integration  loss  for  target  Doppler 
=  [l  +  -N\[?^+N^—  \ 

TT  V  TT  TT 


(A61) 


(A62) 


(A63) 


(A64) 


(A65) 


(A66) 
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(A67a) 


For  DPCA,  the  difference  signals  for  target  (!’)  and  clutter  (C)  are 
Aj-  =  A  —  2B  cos  2d, 

Ac  =  4-2Bexp-is  ^  (ViA/Z))^  (A67b) 

The  results  for  the  two  candidate  systems  discussed  in  Section  3  were  calculated  using 
Eq.  (A60).  To  indicate  the  adequacy  of  this  approximation,  we  can  define  correction 
factors  for  the  target  (T)  and  clutter  (C)  as 

corr  (T)  =  ),  (A68a) 

corr  (C)  =  (A68b) 

The  correction  factors  are  the  ratios  of  the  exact  results,  Eq.  (A30),  for  target  [Eq.  (A36)] 
and  clutter  [Eq.  (A35)],  and  the  approximate  expressions  leading  to  Eq.  (AGO). 

For  the  parameters  discussed  in  Section  3,  the  correction  factors  as  a  function  of  to  are 
shown  in  Figure  15.  It  can  be  seen  that  the  correction  factors  are  essentially  unity  except 
for  small  decorrelation  lengths  at  L-band  where  they  tend  to  further  degrade  performance. 

The  form  of  the  antenna  improvement  factor,  Eq.  (A64),  employed  in  Eq.  (AGO)  is  based 
on  the  simple  estimate  of  Eq.  (A57).  This  is  adequate  if  ^  is  sufficiently  close  to  zero  but 
can  be  seriously  in  error  over  part  of  the  Doppler  spectrum.  As  illustrations  of  this  effect. 
Figure  IG  shows  the  antenna  improvement  factor  as  predicted  by  Eq.  (A52)  and  as  given 
by  Eq.  (AG4)  for  the  three  phase  center  implementation  of  the  two  systems  discussed  in 
Section  3.  It  is  clear  that  both  systems  deviate  significantly  from  the  simple  expression, 
Eq.  (AG4).  Note  that  ^  =  0.072  for  the  L-band  system  and  ^  =  0.064  for  the  UHF  system. 
The  clutter  errors,  Figure  15,  for  large  values  of  £<,  are  due  to  the  inaccurate  calculation 
of  the  antenna  improvement  factor.  The  first  error  was  that  Eq.  (A64)  was  used  with 
V  V  —  Vex  [cf.  Eqs.  (A57)  and  (A58b)]  leading  to  overestimates  of  1.4  dB  at  L-bmid  and 
0.25  dB  at  UHF.  In  addition,  Figure  16  implies  an  additional  0.6  dB  error  at  L-band  and 
about  —0.1  dB  at  UHF.  The  total  errors  are  approximately  2  dB  at  L-band  and  0.15  dB 
at  UHF. 

Finally,  it  should  be  noted  that  the  form  of  Eq.  (AGO)  ignores  the  Doppler  shifts  induced 
by  the  phase  only  contribution  to  the  ionospheric  transfer  function.  As  discussed  in  Ap¬ 
pendix  B,  this  shift  cai.  be  significant.  The  actual  antenna  improvement  factor  that  should 
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Figure  15.  Correction  factors  for  target  (T)  and  clutter  (C)  versus  £q:  (a)  L-band, 
(b)  UHF. 
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be  used  in  Eq.  (AGO)  should  be  evaluated  at  the  value  of  ^  given  in  Eq.  (A59).  Since  the 
term  arising  from  C2  is  a  Gaussian  random  variable,  the  target  can  easily  be  shifted  such 
that  the  interfering  clutter  increases  and  the  antenna  improvement  factor  decreases. 
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APPENDIX  B 

CHARACTERIZATION  OF  THE  ENVIRONMENT 


The  presence  of  an  electron  density  influences  the  propagation  of  an  electromagnetic  wave 
by  introducing  an  additional  (over  vacuum)  differential  phase  change  per  unit  length  of 
the  propagation  path,  given  by 

^  =  reAne(r,<)  (Bl) 

where  (f)  is  the  phase,  rg  is  the  classical  radius  of  the  electron  (rg  =  2.82  x  10“*®  m)  and  Ug 
is  the  local  electron  number  density.  This  assumes  that  the  electron  density  is  not  so  high 
that  absorption  dominates.  The  total  electron  density  can  be  partitioned  into  an  average 
value,  <  Ug  >,  and  a  fluctuating  contribution,  ^ng(r,  t),  as 

ng(r,i)  =  <  rig  >  +  ^ng(r,f).  (B2) 

The  average  value  leads  to  total  electron  content  (TEC)  effects  and  as  can  be  seen  from 
Eq.  Bl,  will  be  frequency  dependent.  For  a  generic  radar  system,  the  TEC  effects  would 
affect  ranging  since  the  coded  pulse  has  some  finite  bandwidth  B.  The  effects  include 
spreading  the  range  resolution  cell  and  introducing  phase  errors  for  the  signal  contribution 
of  individual  scatterers.  However,  the  systems  of  interest  have  reasonably  small  bandwidths 
(megahertz)  and  moderately  high  frequencies  (\'HF  and  above)  so  that  these  effects  only 
become  important  for  severe  environments  where  other  effects  have  already  significantly 
degraded  performance.  For  this  report,  we  will  not  consider  frequency  effects  within  the 
bandwidth  of  the  radar. 

The  effects  of  the  fluctuating  part  can  be  viewed  as  follows  (Refs.  33-34).  Consider  an 
electromagnetic  wave  incident  on  an  extended  region  of  disturbed  electron  density  such  as 
the  ionosphere.  Initially,  the  wave  experiences  phase  changes  that  are  random  from  one 
location  to  another.  These  random  phase  fluctuations  produce  distortions  in  the  wave  front 
which  in  turn  leads  to  diffraction.  As  the  perturbed  wave  propagates  into  the  disturbed 
region,  it  can  break  up  into  additional  ray  paths  that  can  further  diffract  and  scatter. 
Consequently,  if  we  consider  a  point  source  on  one  side  of  a  perturbed  region  and  some 
point  (scatterer)  on  the  other,  the  electromagnetic  field  that  reaches  the  point  may  have 
propagated  along  a  number  of  more  or  less  independent  paths.  The  total  field  present,  being 
a  coherent  sum  of  the  contributions  from  different  paths,  can  have  significant  differences  in 
phase  and  amplitude  from  that  given  by  free  space  propagation.  We  can  define  a  transfer 
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function,  h,  from  a  source  on  one  side  of  the  pertiurbing  media  (r+)  to  a  point  on  the 
other  side  (r_)  as  the  ratio  of  the  actual  electric  field  present  (E)  to  the  free  space  electric 
field  (E,), 

/t(r+,r_,t)  =  E{r+,r-,t)/Eoir+,r-).  (B3) 

Here,  the  time  t  refers  to  changes  in  the  media  which  occur  on  scales  much  slower  than 
the  time  scale  set  by  the  bandwidth  of  the  radar,  1/B. 

The  actual  characterization  of  the  medium,  ne(r,  f),  is  never  known.  Instead,  the  medium 
is  described  in  statistical  terms.  Likewise,  the  solution  to  the  transfer  function  is  not 
available  so  it  is  also  described  statistically.  The  second  order  statistics  are  determined  by 
the  mutual  coherence  function,  F i ,  defined  as 

r,(ri+,ri_,r2+,r2_,ti,t2)  =  (/j(ri+,ri_,fi)/t*(r2+,r2-,t2)).  (B4) 

This  is  the  primary  quantity  o5  interest  in  the  study  of  communication  channels.  However, 
since  radar  employs  two  way  paths  through  the  media,  the  primary  quantity  of  interest  is 
the  fourth  order  statistics,  given,  in  an  obvious  notation,  by 

r2(l,2,3,4)  =  {h{l)h(2)h*{3)h*(4)).  (B5) 

Approximate  analytic  solutions  are  available  for  Fi  in  terms  of  the  statistics  of  6ne(r,  t) 
(Ref.  35).  A  general  result  for  F2  is  not  available.  However,  it  has  been  determined  that 
under  strong  scattering  conditions  (Refs.  36,  37)  the  channel  transfer  function  can  be 
usefully  partitioned  as 

h  =  hs  exp[i<f)i,l  (B6) 

where  (})i  is  a  Gaussian  process  and  the  real  and  in.aginary  parts  of  /i,,  are  zero-mean, 
equal  variance,  uncorrclated  Gaussian  random  variables. 

In  this  formulation,  the  total,  sight-path-integrated  phase  (i.e.,  from  Eq.  Bl)  is  also  parti¬ 
tioned  as  (f>  =  </>,,  -1-  </>/^.  This  is  done  in  such  a  manner  (Ref.  37)  that  (j)[,  is  spatially  slowly 
varying  and  therefore  does  not  cause  significant  scattering  and  multipath  effects.  Then  F 1 
can  also  be  partitioned,  at  least  heuristically,  as 

r,(i,2)  =  r,3(i,2)exp[nVL(i)  -<?/d'2))j,  (B7) 

where  Fj,  is  the  nuitual  coherence  function  arising  from  effects  caused  by  In  the  thin 
phase  .screen  approximation,  F|,  is  giv('n  as 

F,,(1.2)  =  exp(CV,(l,2)-0,  (B8) 
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where  is  the  variance,  and  C^,(l,2)  the  covariance,  of 

Most  important  for  present  purposes  is  the  conclusion  that,  since  h,  may  now  be  regarded 
as  a  complex  Gaussian  process,  with  additional  properties  as  cited  above,  under  these 
strong  scattering  conditions  /i,  is  fully  defined  by  its  second  order  statistics,  or  Fia.  The 
needed  relationship  for  r2  can  then  be  given  as 

r2(l,2,3,4)  =  (expi(^L(l)  +  4>l(2)  -  -  <^i(4))) 

X  [r„(i,3)ri,(2,4)  +  r„(i,4)ri,(2,3)].  (B9) 

A  further  simplification  will  be  obtained  from  the  assumption  that  the  scattering  iono¬ 
spheric  medium  evolves  on  short-to-intermediate  time  scales  (i.e.,  seconds-to-minutes)  pri¬ 
marily  through  an  overall  uniform  drift  relative  to  the  radar  sight  path(s).  With  this 
so-called  “Taylor  frozen-in”  hypothesis,  the  medium  will  be  assumed  to  be  time  indepen¬ 
dent  in  its  rest  frame,  which  is  the  coordinate  system  employed  for  the  calculations  of 
Appendix  A.  In  addition,  Fi,,  Eq.  B8,  will  be  approximated  by  assuming  (a)  structure 
occurs  only  in  the  along-track  direction,  (b)  the  two  sides  of  the  perturbing  media  are 
related  by  the  thin  phase  screen  result,  and  (c)  an  effective  decorrelation  length  at  the 
radar  platform,  Cq,  is  introduced.  The  particular  expression  we  will  employ  is 

F„(l,2)  =  exp-(.Yi  -X2)V^^  (BIO) 

where 

X  =  (Bll) 

Here,  refers  to  the  radar  side  of  the  perturbing  media,  X-  to  the  clutter  and  target  side 
and  R^{R^)  is  the  distance  from  the  center  of  the  perturbing  media  to  the  radar  (clatter). 
We  now  turn  to  a  discussion  of  the  relationship  between  ^nd  the  perturbing  media. 

The  electron  power  spectral  density  (PSl)  is  defined  as 

(6n,(r)^ne(r'))  =  J  (B12) 

Experimental  evidence  suggests  that  a  reasonable  PSD  has  a  power  law  behavior  of  /v 
for  large  K  and  is  syminetiic  around  the  ambient  magnetic  field  direction.  If  w(>  assume  a 
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K  "  behavior,  a  convenient  power  law  PSD  is  (Ref.  33) 

~r~ '  ' 

_ 1 _ 

[I  +  LlKl  +  L\K]  +  LlK-i]' 

X  Jf,  ^/3[1  +  LlKl  +  L\Kl  +  Llh'l]  *  j  (B13) 

where 

^  =  iilLo.  (B14) 

We  have  taken  the  ambient  magnetic  field  direction  in  the  y-direction  to  more  easily 
transition  to  the  1-D  case  employed  in  this  paper.  The  quantities  Lo  and  Ly  are  usually 
referred  to  as  outer  scale  lengths,  with  typically 

Ly  >  io.  (B15) 

The  quantity  is  called  an  inner  scale  length  and  is  such  that 

^  <  1.  (B16) 

With  the  PSD  given  in  Eq.  B13,  we  have 

$«„(K)  ~  ~  -1  <  <~  i  (B17) 

Lq  t, 

In  general,  the  axis  of  symmetry  for  the  electron  fluctuation  PSD  is  not  related  to  the 
platform  motion  and  viewing  geometry.  For  this  more  general  case,  the  main  consequence 
for  the  discussions  of  this  work  would  be  to  introduce  a  general  positive  definite  quadratic 
form  to  incorporate  the  decorrelation  lengths  [cf.  Eq.  BIO]  and  range  and  azimuthal 
effects  would  become  mixed.  We  here  restrict  ourselves  to  the  simpler  1-D  situation  since 
it  contains  all  the  essential  elements  of  the  analysis. 

As  a  ray  traverses  the  disturbed  ionosphere  (here  referred  to  as  the  z-direction),  the  phase 
will  acquire  a  fluctuating  contribution,  (see  Eqs.  B1  and  B2] 

rip  /2 

S(j>{rj_)  =  \re  dz6ne{rx,z)  (B18) 

J  —  ty  /2 


48 


where  Cp  is  the  thickness  of  the  ionosphere.  The  correla  a  function  for  these  phase 
fluctuations  is  then 


2n^  2  f 


g«Kj.-(ri-r2 


)j. 


(B19a) 


or,  for  Ly  — >  oo, 


{S(t>{xi)8<l>{x2))  =  J 


2  f  tKr(xi-Z2) 

•f’Kn^(d)J  27r 


X - ^  1  +  LlKl 

[l  +  LlKl]  *  L 

In  the  above,  the  average  square  of  the  phase  fluctuations  is 

O  /I  »  O  —  O  — n  ■  ■  \  ~  / 


(B19b) 


(tJ  =  \/^A^r^4Lo<TNe\/^ 


Kn.B{(3)' 


Note  that  for  n  =  4  and  /?  small,  we  obtain 


al  =  2\'^rlepLo<TNe 


(B20a) 


(B20b) 


Above,  Eq.  B6,  we  separated  the  transfer  function  into  a  Rayleigh  scintillation  contribution, 
hg,  and  a  phase  only  contribution.  The  separation  is  effected  by  separating  the  scale  sizes 
in  the  phase  PSD, 


=  <tIC(K), 

=  -K^) - (/?[1  +  LlK^]^^), 


(B21a) 

(B21b) 


at  the  Rayleigh  wavevector,  Kfi,  determined  by  the  following  condition  on  the  Rytov 
parameter  (Ref.  37): 

0.1  =  <7^  f  ^^C(K)  sin^  (B22) 

Jo  ^ 


wnere 


R*  = 


R^R^ 

R-\.  "b  R— 


(B23a) 

(B23b) 
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Clearly,  Eq.  B22  sets  a  lower  bound  for  before  Rayleigh  scintillations  occur.  With  the 
typical  situation  of  7  and  being  small  and  for  the  n  =  4  PSD,  cr^  must  satisfy 


2  2  1  2 

(B24) 

A  typical  value  for  Lg  is 

Lg  =  10  km 

(B25) 

so  for  the  candidate  systems  discussed  in  Section  3,  we  find 

<xl  (L-band)  >  1970, 

(B26a) 

al  (UHF)  >  700. 

(B26b) 

(Note  that  Ly,  the  outer  scale  length  in  the  magnetic  field  direction,  is  about  150  km.) 


The  important  parameter  that  cheiracterizes  the  radar  performance  is  the  decorrelation 
length  for  the  Rayleigh  contribution  at  the  radar  platform  location.  The  first  step  is  to 
determine,  from  Eq.  B22,  the  relationship  between  cr^  and  Kr.  Once  Kr  is  determined, 
the  effective  decorrelation  length  satisfies 


_  1 
£l  ~  2 


dK 


TT 


K^C{K). 


(B27) 


Figure  17  shows  the  relationship  between  and  Kr  for  the  two  systems  discussed  in 
Section  3,  again  assuming  the  n  =  4  PSD,  Eq.  B25  for  Lq  and  a  value  for  /9  of  10“^.  The 
decorrelation  length  as  a  function  of  is  determined  from  Eq.  B27.  These  relationships  for 
the  two  systems  are  shown  in  Figure  18.  Notice  that  for  on  the  order  of  unity,  the 

decorrelation  length  changes  very  rapidly  with  £0  — ♦  00  as  >  1-  As  an  example, 

we  will  consider  severe  ambient  spread-F.  In  order  to  estimate  a  set  of  parameters  for 
spread  F,  we  employ  the  electron  density,  meeisured  during  PLUMEX  I  (Ref.  38)  and 
assume  locally 


(in,f  ~ 


Integration  of  the  experimental  data  then  yields 
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We  then  have 


<rl  (L-band)  =368, 
<rl  (UHF)  =2880. 


According  to  Eq.  B26a,  the  L-band  value  is  lower  than  that  required  to  produce  Rayleigh 
scintillations.  The  decorrelation  length  for  the  UHF  system  {cr^jcrM  —  2.03)  is  found  from 
Figure  18  to  be 

4  (UHF)  ~  250  m. 


We  next  consider  the  Doppler  spread  produced  by  the  long  wave-length  components  of  the 
disturbing  media.  With  F,  Eq.  A34,  approximated  as  in  Eq.  A51b,  we  have 

fifC 

Cl  =  2T^V^al  /  —K^C(K).  (B23) 

^  TT 

The  only  important  effect  is  in  Z^,  Eq.  A55.  If  we  define 


N^—=^1+^2, 

TT 

8  - 
^2  = 


(B29a) 

(B29b) 

(B29c) 


Pi  and  P2  characterize  the  Doppler  spread  produced  by  the  Ra2deigh  scintillations  and 
the  long  wavelength  structure,  respectively.  As  illustrations.  Figure  19  shows  Pi  and  P2 
for  the  two  systems  discussed  in  Section  3,  assuming  Lg  =  Vtc.  As  can  be  seen,  the  P2 
contribution  is  either  very  smzdl  compared  to  Pi  or  both  are  reasonably  small  so  that  the 
total  Doppler  spread  is  quite  small. 


Finally,  we  consider  the  Doppler  shift  produced  by  the  long  wavelength  component  of  the 
distiu-bing  media.  This  Doppler  shift,  SVd,  is  given  in  Eq.  A59  as 


SVd  = 


A  C2{C)-C2{T) 


2T  TT 

This  is  a  Gaussian  random  variable  with  a  covari«mce  given  by 


(B30) 


{{SVd?)  = 


2tt^  \R 


2  rdK 


C{K)K^ 


l-cosRT^X 

/t— 


(B31) 
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where 


K  = 


(B34) 


and  X  is  the  separation  between  the  target  and  the  competing  clutter.  As  an  illustration  of 
the  magnitude  of  the  potential  Doppler  shift,  we  will  calculate  the  standard  derivation  when 
the  separation  is  the  azimuthal  extent  of  the  radar  footprint,  X  =  A  A.  The  results  for 
the  two  systems  discussed  in  Section  3  are  shown  in  Figure  20.  Recall  that  the  ambiguous 
Doppler  velocities  (\f2T)  for  the  L-band  systems  are  625  m/sec  for  the  three-phase  center 
approach  and  312  m/sec  for  the  two-phase  cente.''  approach  while  the  corresponding  values 
for  the  UHF  system  are  700  m/sec  and  350  m/sec,  respectively. 


55 


1 


1 

I 


Figure  20.  (tsvd  versus  £o'  (a)  L-band,  (b)  UHF. 


56 


DISTRIBUTION  LIST 

DNA-TR.90-78 


DEPARTMENT  OF  DEFENSE 

ASSISTANT  TO  THE  SECRETARY  OF  DEFENSE 
ATTN.  EXECUTIVE  ASSISTANT 

DEFENSE  ADVANCED  RSCH  PROJ  AGENCY 
ATTN.  CHIEF  SCIENTIST 
ATTN:  GSDRALEWINE 

DEFENSE  COMMUNICATIONS  AGENCY 
ATTN;  C4S/SSM 
ATTN:  SSS 

DEFENSE  COMMUNICATIONS  ENGINEER  CENTER 
ATTN;  CODE  R430  BOEHM 

DEFENSE  INTELLIGENCE  AGENCY 
ATTN:  DB-TPO 
ATTN:  DC-6 
ATTN:  DIR 
ATTN  DT-IB 
ATTN:  RTS-2B 

DEFENSE  NUCLEAR  AGENCY 
ATTN*  NANF 
ATTN:  NASF 
ATTN:  OPNA 
ATTN:  PRPDRYOHO 
3  CYS  ATTN:  RAAE 

ATTN:  RAAEACHESLEY 
ATTN:  RAAE  G  ULLRICH 
ATTN:  RAAE  K  SCHWARTZ 
ATTN:  RAAE  S  BERGGREN 
ATTN  RAEE 
2  CYS  ATTN  TITL 

DEFENSE  NUCLEAR  AGENCY 

ATTN*  FCPRA  LC  R  HEDTKE 
ATTN:  TDNM 

2  CYS  ATTN:  TDTTWSUMMA 

DEFENSE  TECHNICAL  INFORMATION  CENTER 
2  CYS  ATTN.  DTIC/FDAB 

JOINT  DATA  SYSTEM  SUPPORT  CTR 
ATTN:  JNSV 

NATIONAL  SECURITY  AGENCY 
ATTN  LPLUSWICK 

STRATEGIC  AND  THEATER  NUCLEAR  FORCES 
ATTN  DRESEVIN 
ATTN  DRSCHNEITER 
ATTN.  GWINNEBURG 

STRATEGIC  DEFENSE  INITIATIVE  ORGANIZATION 
ATTN  EN 

ATTN  EN  LTCC  JOHNSON 
ATTN*  PTPLSTOESSELL 
ATTN  TN  DR  M  GRIFFIN 


THE  JOINT  STAFF 

ATTN:  JKC  (ATTN*  DNA  REP) 

ATTN:  JKCS 

ATTN:  JLWT  THREAT  ANALYSIS 
ATTN:  JPEM 

THE  JOINT  STAFF 
ATTN:  J6 

U  S  NUCLEAR  CMD  &  CENTRAL  SYS  SUP  STAFF 
ATTN:  SABHSEQUINE 

DEPARTMENT  OF  THE  ARMY 

ARMY  LOGISTICS  MANAGEMENT  CTR 
ATTN:  DLSIE 

HARRY  DIAMOND  LABORATORIES 
ATTN.  SLCIS-IM-TL 

U  S  ARMY  ATMOSPHERIC  SCIENCES  LAB 
ATTN:  SLCAS-AE  DR  NILES 
ATTN;  SLCAS-AR  DR  H  HOLT 

U  S  ARMY  COMMUNICATIONS  R&D  COMMAND 
ATTN;  AMSEL-RD-ESA 

U  S  ARMY  ENGINEER  DIV  HUNTSVILLE 
ATTN;  PRESTON  J  KISS 

U  S  ARMY  FOREIGN  SCIENCE  &  TECH  CT'? 

ATTN:  AIFRTA 

U  S  ARMY  MISSILE  COMMAND/AMSMI-RD-CS-R 
ATTN;  AMSMI-RD-CS-R 

U  S  ARMY  NUCLEAR  &  CHEMICAL  AGENCY 
ATTN.  MONA-NUDBASH 

U  S  ARMY  NUCLEAR  EFFECTS  LABORATORY 
ATTN:  ATAA-PL 
ATTN.  A1AA-TDC 

ATTN:  ATRC-WCC  LUIS  DOMINGUEZ 

U  S  ARMY  STRATEGIC  DEFENSE  CMD 
ATTN:  CSSD-H-LSBCARRUTH 
ATTN*  CSSD-H-SA 
ATTN:  CSSD-H-SAV 
ATTN;  CSSD-SA-EV  RON  SMITH 

U  S  ARMY  STRATEGIC  DEFENSE  COMMAND 
ATTN.  CSSD-GR-SW  DICKINSON 

USA  SURVIVABILITY  MANAGMENT  OFFICE 
ATTN  SLCSM-SE  J  BRAND 

DEPARTMENT  OF  THE  NAVY 

COMMAND  &  CONTROL  PROGRAMS 
ATTN:  OP  941 

DEPARTMENT  OF  THE  NAVY 
ATTN  JCMG-707 


Dist-1 


DNA-TR-90-78  (DL  CONTINUED) 


NAVAL  AIR  SYSTEMS  COMMAND 
ATTN  PMA271 

NAVAL  ELECTRONICS  ENGRG  ACTVY  PACIFIC 
ATTN  CODE  250 

NAVAL  RESEARCH  LABORATORY 

ATTN.  CODE  2000  J  BROWN 
ATTN.  CODE  2627 

ATTN:  CODE 4121.8 HHECKATHORN 
ATTN-  CODE  4183 
ATTN  CODE  4701 
ATTN:  CODE  4720  J  DAVIS 
ATTN-  CODE  4780  BRIPSN 
ATTN:  CODE  4780  DR  P  BERNHARDT 
ATTN  CODE  4780  J  HUBA 
ATTN  CODE  4785  P  RODRIGUEZ 
ATTN-  CODE  5300 
ATTN-  CODE  5326  6  A  ANDREWS 
ATTN.  CODE  5340  E  MOKOLE 
ATTN:  CODE  8344  M  KAPLAN 
ATTN-  JACOB  GRUN 

NAVAL  SURFACF  WARFARE  CENTER 
ATTN-  CODEH-21 

NAVAL  TECHNICAL  INTELLIGENCE  CTR 
ATTN-  DA44 

NAVAL  UNDERWATER  SySTEMS  CENTER 
ATTN  CODE  341 IJKATAN 

OFFICE  OF  CHIEF  OF  NAVAL  OPERATIONS 
ATTN.  OP  654 
ATTN-  OP  9410 

SPACE  &  NAVAL  WARFARE  SYSTEMS  CMD 
ATTN-  CODE  3101  T  HUGHES 
ATTN  PD50TD 
ATTN.  PD50TD1  G  BRUNHART 
ATTN  PME  106*4  S  KEARNEY 
ATTN  PME-106F  WDIEDERICH 

THEATER  NUCLEAR  VMRFARE  PROGRAM  OFC 
ATTN.  PMS-42331FO  SMITH 

DEPARTMENT  OF  THE  AIR  FORCE 

AFIA/INKS 

ATTN  AFIA/INKS  MAJ  SCHROCK 
ATTN  AFIA/INKS 

AIR  FORCE  CTR  FOR  STUDIES  &  ANALYSIS 
ATTN  AFCSA/SASC 

AIR  FORCE  ELECTRONIC  WARFARE  CENTER 
flTTM  LJ  M  MONEELY 
ATTN  SAZ 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RSCH 
ATTN  AFOSR/NP 


AIR  FORCE  PHILLIPS  LABORATORY 
ATTN:  J  KLOUBACHAR 
ATTN  OP/WBLUMBERG 
ATTN.  SANTI  BASU 

AIR  FORCE  SYSTEMS  COMMAND 
ATTf'-  XTTW 

AIR  FORCE  SYSTEMS  COMMAND 
ATTN.  J  COLYE' 

AIR  UNIVERSITY  LIBRARY 
ATTN  AUL-LSE 

HQ  AWS,  DET  3  (CSTC/WE) 

ATTN  WE 

NATIONAL  TEST  FACILITY 

ATTN  NTB/JPO  DR  C  GIESE 

PHILLIPS  LABORATORY 
ATTN-  NTCA 

ATTN  NTCTS  LTC  C  ROYER 
ATTN  NTN 

STRATEGIC  AIR  COMMAND/XRFS 
ATTN  IN 
ATTN.  XRFS 

DEPARTMENT  OF  ENERGY 

EG&G,  INC 

ATI  IN.  D  WRIGHT 

LAWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN-  L-97TDONICH 

SANDIA  NATIONAL  LABORATORIES 
ATTN-  P  L  MATTERN  8300 

SANDIA  NATIONAL  LABORATORIES 

ATTN.  A  D  THORNBROUGH  400 
ATTN  CODE  9014  R  BACKSTROM 
ATTN  D  DAHLGREN  6410 
ATTN  DIV  2344  ROBERT  M  AXLINE 
ATTN  ORG  9110  G  CABLE 
ATTN  ORG  9110  WO  BROWN 
ATTN  TECH  LIB  3141 

OTHER  GOVERNMENT 

CENTRAL  IN  FELLIGENCE  AGENCY 
ATTN.  OSWR/NED 
ATTN  OSWR/SSD  FOR  L  BERG 

DEPARTMENT  OF  COMMERCE 
ATTN-  G  REEVE 
attm  I  unccft^cvco 

I  .^1  4\/t  t  IflU  f 

ATTN  WUTLAUT 

U  S  DEPARTMENT  OF  STATE 
ATTN  PM/TMP 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORP 

ATTN  A  MORSE 
ATTN.  BRIAN  PURCELL 
ATTN:  C  CREWS  92957 
ATTN:  CRICE 
ATTN  OR  JM  STRAUS 
ATTN.  OUGHT 
ATTN:  IGARFUNKEL 
ATTN-  JKLUCK 
ATTN.  J THACKER 
ATTN  MROLENZ 

AT&T  BELL  LABORATORIES 

ATTN:  DENIS  SLONGO 
ATTN:  JOSEPH  A  SCHOLL 
ATTN-  N  BEAUCHAMP 

ATLANTIC  RESEARCH  SERVICES  CORP 
ATTN:  R  MCMILLAN 


EOS  TECHNOLOGIES  JNC 
ATTN.  B  GABBARD 
ATTN:  RLELEVIER 

FORD  AEROSPACE  CORPORATION 
ATTN:  PATRICIA  BIRDWELL 

GENERAL  ELECTRIC  COMPANY 

ATTN  JOSEPH  E  STROSSER 

GENERAL  RESEARCH  CORP  INC 
ATTN-  JEOLL 

GRUMMAN  AEROSPACE  CORP 
ATTN-  JDIGLIO 

HARRIS  CORPORATION 

ATTN:  LYMUEL  MCRAE 

HSS,  INC 

ATTN.  D HANSEN 


ATMOSPHERIC  \hO  ENVIRONMENTAL  RESEARCH  INC 
ATTN  MKO 


INFORMATION  SCIENCE.  INC 
ATTN:  WDUDZIAK 


AUSTIN  RESEARCH  ASSOCIATES 
ATTN.  J  THOMPSON 

AUTOMETRIC.  INC 

ATTN.  C  LUCAS 

BDM  INTERNATIONAL  INC 

ATTN.  W  LARRY  JOHNSON 

BERKELEY  RSCH  ASSOCIATES.  INC 
ATTN.  J  WORKMAN 
ATTN:  NTGLADD 
ATTN:  S  BRECHT 

BOEING  CO 

ATTN:  GHALL 

CALIFORNIA  RESEARCH  &  TECHNOLOGY.  INC 
ATTN.  M  ROSENBLATT 

CHARLES  STARK  DRAPER  LAB.  INC 
ATTN.  ATETEWSKl 

COMMUNICATIONS  SATELLITE  CORP 
ATTN  RICHARD  A  ARNDT 

CONTEL  FEDERAL  SYSTEMS  INC 

ATTN.  CHARLES  BENNINGTON 

CORNELL  UNIVERSITY 

ATTN-  D  FARLEY  JR 
ATTN  M  KELLY 

D^NETICS.  INC 

ATT  T  WILLIAM  D  TEPPER 

ELECTROSPACE  SYSTEMS.  INC 

ATTN-  LINDA  CALDWELL 
ATiN  P  PHILLIPS 


INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN.  E  BAUER 
ATTN.  HWOLFHARD 

JAYCOR 

ATTN-  AGLASSMAN 
ATTN:  J  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
ATTN.  CMENG 
ATTN:  HGTORNATORE 
ATTN-  J  D  PHILLIPS 
ATTN  R  STOKES 

KAMAN  SCIENCES  CORP 
ATTN:  DASiAC 
ATTN.  E  CONRAD 
ATTN:  GDITTBERNER 

KAMAN  SCIENCES  CORPORATION 
ATTN:  BGAMBILL 
ATTN-  DASIAC 
ATTN  R  RUTHERFORD 

LOCKHEED  MISSILES  &  SPACE  CO.  INC 
ATTN:  J  HENLEY 
ATTN-  J  KUMER 
ATTN  R SEARS 

LOCKHEED  MISSILES  &  SPACE  CO.  INC 
ATTN.  CARLCRABILL 
ATTN-  DTRAMPTON 
ATTN.  EMDIMZCELI 

LTV  AEROSPACE  &  DEFENSE  COMPANY 
2  CYS  ATTN.  LIBRARY  EM-()8 

M  I  T  LINCOLN  LAB 

ATTN-  D  TOWLE  L.230 
ATTN.  IKUPIECL  IOO 
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MARTIN  MARIETTA  DENVER  AEROSPACE 
ATTN:  J  BENNETT 
ATTN-  H  VON  STRUVE  111 


MAXIM  TECHNOLOGIESJNC 


ATTN: 

B  PHILLIPS 

ATTN: 

N  CIANOS 

MAXWELL  LABS,  INC 

ATTN: 

KWARE 

MCDONNELL  DOUGLAS  CORPORATION 

ATTN: 

J  GROSSMAN 

ATTN: 

R  HALPRIN 

METATECH  CORPORATION 

ATTN: 

W  RADASKY 

METEOR  COMMUNICATIONS  CORP 

ATTN: 

R  LEADER 

MISSION  RESEARCH  CORP 

ATTN: 

R  ARMSTRONG 

ATTN: 

W  WHITE 

MISSION  RESEARCH  CORP 

ATTN: 

R  L  BOGUSCH 

MISSION  RESEARCH  CORP 

ATTN: 

DAVEGUICE 

MISSION  RESEARCH  CORP 

ATTN- 

BR  MILNER 

ATTN: 

CLONGMIRE 

ATTN: 

D  KNEPP 

ATTN: 

D  LANDMAN 

ATTN: 

F  FAJEN 

ATTN: 

FGUIGLIANO 

ATTN. 

G  MCCARTOR 

ATTN: 

KCOSNER 

ATTN: 

M  FIRESTONE 

ATTN. 

R  BIGONI 

ATTN: 

RDANA 

ATTN: 

R  HENDRICK 

ATTN: 

RKILB 

ATTN: 

SGUTSCHE 

ATTN: 

TECH  LIBRARY 

MITRE  CORPORATION 

ATTN.  M  HORROCKS 
ATTN-  RCPESCiW-34 
ATTN:  W  FOSTER 


MITRE  CORPORATION 

ATTN:  GCAMPARETTO 

NORTHWEST  RESEARCH  ASSOC.  INC 
ATTN:  EFREMOUW 

PACIFIC-SIERRA  RESEARCH  CORP 
ATTN:  E  FIELD  JR 
ATTN:  HBRODE 


PHOTOMETRICS.  INC 

ATTN:  1  L  KOFSKY 

PHOTON  RESEARCH  ASSOCIATES 
ATTN-  DBURWELL 

PHYSICAL  RESEARCH  INC 
ATTN:  WSHIH 

PHYSICAL  RESEARCH  INC 
ATTN-  ACECERE 

PHYSICAL  RESEARCH  INC 
ATTN:  HFITZ 

PHYSICAL  RESEARCH.  INC 
ATTN:  RDELIBERIS 
ATTN:  T  STEPHENS 

PHYSICAL  RESEARCH,  INC 
ATTN:  J  DEVORE 
ATTN:  J  THOMPSON 
ATTN.  WSCHLUETER 

PHYSICS  INTERNATIONAL  CO 
ATTN:  C  GILMAN 

R  &  D  ASSOCIATES 
2  CYS  ATTN:  B  LAMB 
ATTN:  GHOYT 
5  CYS  ATTN:  L  DEREAD 

R  &  D  ASSOCIATES 

ATTN;  D  CARLSON 
ATTN  J  THOMPSON 

RAND  CORP 

ATTN:  C  CRAIN 
ATTN:  EBEDROSIAN 

RAND  CORP 

ATTN:  B  BENNETT 

RJO  ENTERPRISES/POET  FAC 

ATTN.  STEVEN  KRAMER 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  R POTTER 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN-  C  SMITH 
ATTN:  D  SACHS 
ATTN.  LLINSON 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN-  HSUNKENBERG 
ATTN:  LIBRARY 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN-  SROSENCWEIG 

SPARTA  INC 

ATTN:  DDEAN 
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SRI  INTERNATIONAL 

ATTN:  R  LIVINGSTON 
ATTN:  RTTSUNODA 
ATTN-  WCHESNUT 
ATTN:  WJAYE 

STEWART  RADIANCE  LABORATORY 
ATTN:  RHUPPI 

TELECOMMUNICATION  SCIENCE  ASSOCIATES 
ATTN:  R  BUCKNER 

TELEOYNE  BROWN  ENGINEERING 
ATTN:  J  FORD 
ATTN-  J  WOLFSBERGER  JR 
ATTN:  N  PASSING 

TOYON  RESEARCH  CORP 
ATTN:  J  ISE 

TRW  INC 

ATTN:  ED  SIMMONS 


TRW  SPACE  &  DEFENSE  SECTOR  SPACE  & 
ATTN:  DM  LAYTON 

USER  SYSTEMS,  INC 

ATTN:  SWMCCANDLESSJR 

UTAH  STATE  UNIVERSITY 
ATTN:  K  BAKER 
ATTN:  L JENSEN 

VISIDYNE,  INC 

ATTN:  J  CARPENTER 

FOREIGN 

FOA2 

ATTN:  BSJOHOLM 

DIRECTORY  OF  OTHER 

BOSTON  UNIVERSITY 

ATTN:  MICHAEL  MENDILLO 
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